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Introduction 


Over  the  past  several  years  S-CUBED,  under  DARPA  support,  has 
been  developing  a  comprehensive  new  seismic  yield  estimation  system  for 
underground  nuclear  explosions.  The  principal  objective  of  this  program 
has  been  to  implement  a  flexible,  interactive  software  system  in  which 
yield  estimates  based  on  a  wide  variety  of  different  seismic  magnitude 
measurements  can  be  efficiently  determined,  merged  with  all  available 
information  regarding  the  test  location  under  consideration  and  statistically 
combined  to  obtain  both  a  unified  seismic  estimate  of  explosion  yield  and 
quantitative  measures  of  the  uncertainty  in  that  estimate.  A  preliminary 
prototype  version  of  a  system  designed  to  achieve  the  above  objective,  des¬ 
ignated  the  Yield  Estimation  System  (YES),  has  been  implemented  in  a 
Sun  color  workstation  (SPARCStation)  environment  at  the  DARPA  Center 
for  Seismic  Studies  (CSS).  This  document  provides  a  user’s  guide  which  is 
applicable  to  the  current  version  of  that  system.  A  comprehensive  overview 
of  the  system  software  design  and  technical  specifications  is  provided  in  a 
recent  report  by  Murphy  et  al.  (1991). 

In  its  current  configuration,  the  system  is  applicable  to  underground 
nuclear  explosions  at  the  Soviet  Shagan  River  and  Novaya  Zemlya  test 
sites  and  encompasses  a  database  of  more  than  10,000  digital  seismograms 
recorded  at  stations  of  the  USAEDS,  GDSN,  CDSN  and  IRIS  networks 
from  explosions  at  these  two  sites.  For  both  test  areas,  information  regard¬ 
ing  the  explosion  .source  environment  is  presented  to  the  analyst  in  the  con¬ 
text  of  SPOT^^  satellite  images  of  the  sites,  together  with  associated 
surface  and  subsurface  geologic  information  and  DMA  topographic  data. 
The  on-line  database  for  YES  also  contains  a  wide  variety  of  tabular  infor¬ 
mation,  including  complete  event  and  station  location  files  containing  both 
classified  and  unclassified  locations,  standard  travel-time  tables  for  the 
>mr oiiant  phases  (Herrin  et  al.,  1968),  propagation  path  and  station  correc- 
'h'ivs  lor  use  in  magnitude  determinations  and  a  comprehensive  instrument 
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response  database.  The  instrument  database  is  a  compilation  ol  information 
collected  from  the  USGS,  AFTAC,  IRIS  and  others  which  details  char¬ 
acteristics  of  the  recording  instruments  as  a  function  of  station,  channel  and 
date.  There  are  currently  750  different  response  histories  in  this  section  of 
the  database. 

The  user  interface  to  the  system  has  been  designed  to  be  completely 
menu-driven  and  mouse-activated  and  requires  no  keyboard  entry  by  the 
operator  once  the  system  has  been  acces.sed.  This  graphical  user  interface 
has  been  built  using  the  X  Window  system.  X  was  designed  specifically  to 
allow  hardware  independence,  to  foster  ease  in  porting  applications  to 
machines  other  than  those  for  which  they  were  developed,  and  to  permit 
running  of  applications  on  one  computer  while  displaying  their  output  on 
another,  even  if  the  computers  are  of  different  manufacture.  YES  was 
designed  and  implemented  on  Sun  SPARC  computers  (e.g.,  SPARCSta- 
tions),  but  transfer  to  any  other  system  that  supports  X  and  UNIX  could  be 
accomplished  easily.  The  X  interface  has  been  written  using  the  X  toolkit. 
The  X  toolkit  enforces  an  object-oriented  approach  to  programming  by 
combining  the  windows  and  the  operations  on  the  windows  into  “widgets." 
YES  uses  several  widget  sets  including  the  Motif  widget  set  from  the  Open 
Software  Foundation,  the  X  widget  set  from  Hewlett  Packard,  the  Athena 
widget  set  from  MIT,  graphics  widgets  written  by  Teledyne-Geotech  and 
SAIC,  and  special  purpose  widgets  for  YES  written  by  S-CUBED. 

YES  con.sists  of  a  hierarchy  of  programs.  At  the  lop  level  is  a  mas¬ 
ter  program  whose  primary  function  is  to  start  the  other  anlaysis  programs. 
The  top  level  program  allows  the  user  to  select  the  lest  site,  event  and  seis¬ 
mic  phase  to  be  processed,  as  well  as  the  analysis  tools  to  apply  to  the  seis¬ 
mic  data  corresponding  to  that  selection.  It  then  starts  the  appropriate 
analysis  programs  with  the  proper  initialization  values.  All  the  required 
software  has  been  written  in  either  the  C  or  FORTRAN  programming  lan¬ 
guages.  C  was  chosen  because  programmer’s  calls  to  the  X  Window  Sys¬ 
tem  procedures  are  in  C,  and  becau.se  C  is  well-suited  to  the  design  of 
complex  sy.sicms  with  a  variety  of  data  structures.  FORTRAN  was  u.sed 
because  it  is  optimal  for  certain  types  of  computational  analysis,  and 
because  its  use  permitted  the  inclusion  of  many  previously  existing  analy- 
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sis  tools  iiiK'  y’ES.  As  a  result,  the  top  level  programs,  interactive  modules, 
datahas',-  ir  ■'!■,,  e,  and  graphics  routines  are  written  in  C,  and,  when  appropri¬ 
ate.  these  n  utines  call  FORTRAN  subroutines. 

This  document  provides  a  user’s  guide  and  reference  manual  for  the 
operation  of  the  current  version  of  the  YES  system.  In  Section  II,  the  capabili¬ 
ties  and  functionality  of  the  .system  are  graphically  illustrated  through  displays 
of  the  screens  encountered  by  an  analyst  in  a  typical  processing  session  for  a 
selected  explosion  (i.e.,  the  Soviet  JVE  explosion  of  14  September  1988).  This 
is  followed  in  Section  III  by  a  brief  summary  of  the  current  status  of  the  YES 
system,  .A  detailed  script  of  the  operator  actions  required  to  reproduce  this 
sample  session  is  provided  in  Appendix  B.  This  is  preceded  in  Appendix  A  by 
an  overview  of  the  software  and  hardware  required  to  implement  YES,  as  well 
as  detailed  startup  instructions  for  the  system.  Appendix  C  provides  a  complete 
description  of  all  the  currently  available  system  options,  including  information 
regarding  their  access  and  application  from  the  menus  contained  in  the  individ¬ 
ual  analysis  modules.  This  is  followed  in  Appendix  D  by  a  top  level  schematic 
diagram  of  the  main  modules  of  the  system  and  their  interdependencies. 
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Overview  of  Selected  VES 
Capabilities 


In  this  section,  some  of  the  capabilities  and  functionality  of  the 
YES  system  will  be  graphically  illustrated  through  displays  of  the  screens 
encountered  by  an  analyst  in  a  typical  processing  session  for  a  selected 
explosion.  For  this  example,  the  unclassified  data  from  the  Soviet  JVE 
explosion  of  14  September  1988  will  be  analyzed  and  used  to  illustrate  var¬ 
ious  features  of  the  system.  A  detailed  script  of  the  operator  actions 
required  to  reproduce  this  sample  session  is  provided  in  Appendix  B.  This 
is  followed  in  Appendix  C  by  a  complete  description  of  all  the  system 
options  which  are  available  in  the  current  version. 


Main 

Menu 


Figure  1. 

Main  menu 
structure 
for  YES. 


As  has  been  noted  previously,  a  distinguishing  characteristic  of  the 
YES  is  that  it  is  completely  menu-driven  and  mouse-activated  and  requires 
no  keyboard  entry  by  the  user.  The  lop  level  menu  providing  access  to  the 
system  is  shown  in  Figure  1 .  It  consists  of  six  buttons  which  can  be  used  to 
initiate  (SITE,  EVENT,  PHASE,  FUNCTION)  or  tenninaie  (QUIT)  action 
w'ithin  the  system  or  to  view  on-line  information  regarding  the  operating 
characteristics  and  parameters  of  the  system  (HELP).  Choosing  any  of 
these  buttons  causes  a  series  of  pulldown  menus  to  be  activated  as  illus¬ 
trated  in  Figure  2.  In  this  example,  the  SITE  button  has  been  activated  to 
select  the  Soviet  Shagan  River  (Balapan)  lest  site,  and  the  EVENT  button 
has  been  activated  to  .select  the  JVE  explosion  which  was  detonated  at  that 
test  site  on  14  September  1988.  Once  the  test  site  and  event  have  been 
selected,  the  remaining  analyst  interaction  with  the  system  is  initiated 
through  the  PHASE  and  FUNCTION  buttons,  llie  PHASE  button  provides 
the  analyst  with  the  capability  to  choose  from  among  the  six  different  seis¬ 
mic  phases  listed  in  Figure  2  for  which  digital  wavelorm  data  are  currently 
available  on  the  system.  The  FUNCTION  button  provides  access  to  the 


4 


t  h'eniew  oi  Sele .  ed  I^ES  Capui  ilities 

seven  pri.ncipai  computational  and  analysis  modules  which  permit  the  ana- 

\y<{  !■' 


view  the  seismic  data  within  the  context  of  the  available  informa¬ 
tion  regarding  the  specific  test  location  under  investigation  Satellite 
Imaf’e,  World  Map) 

interact  with  the  recorded  seismic  data  to  process  it  and  extract  the 
various  magnitude  measures  of  interest  (Analyst  Station,  Magnitude 
Measurement) 


formally  combine  the  seismic  measures  of  source  size  to  obtain  an 
optimum  measure  of  explosion  yield  and  quantitative  measures  of 
the  uncertainty  in  that  estimate  (Yield  Estimation) 


Figure  2. 

Illustration  of 
pulldown  access 
to  individual 


SitelllEvent  Phase 


Fb2»e;  not  set 


menus. 


Events  (1988) 


February  13 


Aprils 


May  4 


FuncUon  Selection 


Satellite  Image 
World  Map 
Analyst  Station 
Magnitude  Measurement  - 
Yield  Estimation 
Statistical  &unmary 
Spreadsheet 
Event  Summa 


[lecember  17 


Statistically  assess  the  results  with  respect  to  any  existing  treaty 
thresholds  or  other  yield  levels  of  particular  interest  (Statistical 
Summary,  Spreadsheet) 


view  a  one  page  summary  of  the  results  obtained  by  applying  the 
system  to  the  data  recorded  from  the  sele^-t-d  explosions  (Event 
Summary). 
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Satellite 

Image 


Figure  3. 

SPOT  satellite 
image  of  the 
Shagan  River 
test  site  with 
superimposed 
locations  of 
the  historical 
explosions 
(squares)  and 
current  event 
(diamond). 


Having  specified  a  test  sue  anc  a  paiiiciiiar  expiosi  .1  (in  this  case 
the  JVE  event),  a  lypicai  anal-  sis  m  qu;  n  a-  ■*  ri.M  i'  'l"  1  svitn  the  ,.*ction 
of  the  Satellite  Image  option  from  the  ^  t  '  e  Us  n.  rni,  v-in>  h  >  nig-;  to 
the  screen  the  SPOT  satellite  image  display  of  the  Sliagan  River  test  site 
shown  in  Figure  3.  In  this  initial  display  of  the  test  site  information  inter¬ 
face,  the  locations  of  previous  explosions  a!  tl.is  •  iu.  so-  shown  us  color- 
coded  square  overlays,  with  the  current  event  highlighted  by  a  yellow  dia¬ 
mond.  In  this  case,  the  different  colors  are  used  to  differentiate  those  explo¬ 
sions  about  which  the  Soviets  have  published  data  in  the  open  literature 
(blue)  from  those  for  which  only  seismic  information  is  available  (red).  It  is 
important  to  note  that  this  is  not  merely  a  static  display,  but  that  in  fact  the 
image  and  the  overlays  to  it  are  formally  tied  to  an  extensive  on-line  data¬ 
base  of  supplementary  information.  Thus,  for  example,  in  this  figure  the 
operator  has  clicked  on  one  of  the  blue  squares  (arrow),  which  has  initiated 
a  process  by  which  available  information  about  that  event  has  been 
extracted  from  the  database  and  displayed  on  the  information  line  below 
the  menu  buttons,  indicating  that  this  explosion  was  detonated  at  Shagan 
on  10  December  1972  at  the  .specified  latitude  and  longitude  and  that 
Bocharov  et  al.  (1989)  have  reported  the  depth  as  478  m  and  the  yield  as 
140  kt. 


Function ;  Help  |[Qi^ 


Images!  iLocationsjljOvcrlays 
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Sham  10  Deocaibcr  1972:  S0.03%N  79.0IO(tE  depth  -  478.0  m 
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Full  res  Fur  purposes  ol  display,  the  satellite  image  shown  in  Figure  3  has 

l  i  ’  ssod  to  an  eiTective  resirlulion  of  about  30  m  so  that  the  entire 
test  site  can  be  viewed  on  a  single  screen.  However,  the  data  corresponding 
to  the  full  10  in  resolution  of  the  SPOT  image  can  also  be  viewed  by  acti¬ 
vating  the  appropriate  button  (Full  res)  on  the  right  hand  margin  of  this  dis- 
plav  and  simply  pointing  with  the  mouse  to  any  location  on  the  image.  For 
example.  Figure  4  shows  the  full  resolution  sub-image  corresponding  to  the 
location  of  the  cratering  explosion  of  15  January  1965  which  dammed  the 
Shagan  River,  producing  the  prominent  lake  in  the  southeast  quadrant  of 
Figure  3  Another  feature  provided  by  this  image  display  module  is  the 
capability  t.  interactively  adjust  the  contrast  and  brightness  using  the  slider 
bars  located  at  the  bottom  of  the  right  hand  menu  margin.  Figure  5  illus¬ 
trates  this  feature  by  way  of  a  comparison  of  the  nominal  display  (right) 
with  that  resulting  from  interactively  reversing  the  contrast  and  decreasing 
the  brightness  (left)  by  repositioning  these  sliders  with  the  mouse.  It  can  be 
seen  that  different  features  are  emphasized  in  these  two  displays,  which 
facilitates  analyst  identification  of  both  geologic  and  man  made  features. 

Surface  Geology  in  addition  to  the  event  locations,  a  number  of  geologic  and  topo- 

and  Topography  graphic  databases  have  been  registered  to  this  SPOT  image  and  are  avail¬ 
able  for  display  from  the  OVERLAYS  and  IMAGES  menus.  For  example, 
Figure  6  shows  the  overlay  of  Leith’s  (1989)  color-coded,  surface  geologic 


Figure  4. 

Full  resolution 
SPOT  satellite 
image  of  the 
region  sur¬ 
rounding  the 
Shagan  River 
cratering 
explosion  of  15 
January  1965. 
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Figure  6. 

SPOT satellite 
image  of  the 
Shagan  River 
test  site  with 
superimposed 
surface  geo¬ 
logic  map 
and  current 
event  loca¬ 
tion. 
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Figure  7. 

Color-coded 
representation 
of  DMA  topo¬ 
graphic  data 
for  the  Shagan 
River  test  site 
with  superim¬ 
posed  topo¬ 
graphic 
contours  and 
current  event 
location. 


I  ■  I  '-i  ,-d  YI'S  Cdp  :l  ilifiev  \ 
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ril'ip  <  !(»■  .1!'.-;,.  where  :tte  ir^cjiiir^r  i  ,  ;c  ..jaeiit  evenr  (square)  is  shown 
in  ■  a  of  the  e'  j  o  \]  '(  riic  units  and  the  trace  of  the 

proininciit  Thinrau  fault  which  intersects  this  portion  of  the  test  site.  The 
Defense  Mapping  Agency’s  topographic  database  for  the  area  is  displayed 
in  color-coded  image  form  in  Figure  7  together  with  overlays  of  the  corre- 
spiv.iding  K'pography  contours  and  current  event  location  selected  from  the 
OVERLAYS  menu. 

The  variation  of  surface  topography,  as  well  as  subsurface  geology 
to  a  depth  of  about  I  km,  along  any  specified  line  across  the  image  can  be 
acce,s,sec!  leing  the  cross-section  button  (Cross  sec)  on  the  right  hand  mar¬ 
gin  of  the  image  display.  This  feature  is  illustrated  in  Figure  8  where  the 
left  hand  panel  shows  the  line  selected  by  the  analyst  by  pointing  with  the 
mouse  to  two  arbitrary  points  on  the  SPOT  image  and  the  right  hand  panel 
shows  the  resulting  vertical  section  through  that  line  which  was  automati¬ 
cally  produced  by  the  system,  together  with  the  locations  and  approximate 
depths  of  penetration  of  explosion  emplacement  holes  encountered  along 
that  line. 
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Figure  8. 

Vertical  subsurface 
section  through  the 
current  event  shot- 
point  along  the  inter¬ 
actively  selected  line 
shown  on  the  SPOT 
image  insert. 


Depth  to  As  with  the  suifacc  typography,  the  variation  in  the  depth  to  any 

Granite  selected  geologic  iia^riacc  can  al,^o  he  exhibiu.d  in  linage  form,  as  Ulus- 

trated  for  the  granite  siii  facc  in  Figure  9  vvheic,  once  again,  the  correspond¬ 
ing  depth  contour,',  .i.i'l  '  'ca';  ■’  iitr  'Vo  w-'cn  overlaid  for 

reference  purposes.  It  is  lir.m, uiaieiy  cvidciu  imm  uii.s  presentation  that. 


Figure  9. 

Color-coded 
representation 
of  depth  to  the 
top  of  the  gran¬ 
ite  surface 
beneath  the 
Shagan  River 
test  site  with 
superimposed 
depth  contours 
and  current 
event  location. 


tmin 


•'  'nr.es 


.  .  V-  ,.!  this  N.iC  hav-. 

'■-■-i '!  •:  .A'.  '  ‘.i  '  initiCL' throiijjhoiil  much  ol' 

ill's  K  .  "1  .  up.iioiiniy  hclow  typical  cspli's;.  u  cmpiacemcnt  depths. 


World  liavitiL  completed  an  iniiia!  u  view  o!  the  available  information 

rcf;  "i.iiy  i.iC  explosion  test  environmept,  the  analyst  can  next  proceed  to 
an  exairiination  of  the  ciare.spondine  recorded  seismic  data.  As  an  initial 
s'l  1  ,  p  -i  ."cw  oi  the  locadons  i  f  stations  lor  which  data  from  the 
sol'  cted  (  xploMon  are  available  iip  the  .systeiii  can  be  obtained  by  selecting 
the  >>■/’(/'/  Map  option  iVcim  the  FLiNCi  U  )NS  menu.  The  resulting  di.splays 
I  1  •  ■  ■  s.  siiiic  1  and  regii  nai  I  „  jii'.ases  are  shown  in  Figure  10, 

plotted  on  a/umithal  equidistant  projections  of  the  globe  (A  <  100°)  cen¬ 
tered  on  the  Shagan  River  lest  site  The  .STATION  menu  item  within  this 
modtilc  provides  the  user  with  the  capability  to  select  any  one  of  these 
recording  stations  from  a  li.st  and  subsequently  cau.se  a  straight  line  to  be 
drawn  on  the  screen  between  the  .source  and  the  receiver  shown  in  the  fig¬ 
ure.  In  this  iransformaiion.  the  straight  line  coiresponds  to  the  great  circle 
path  and  thus  provides  the  analyst  with  a  view  of  the  surface  projection  of 

Figure  10. 

World  Map  projections  (A  <  100°)  showing  locatiims  of  stations  for  which  digital 
teleseismic  P  (left)  and  ( right)  data  are  available  for  the  current  event. 
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Figure  11. 

Analyst  sta¬ 
tion  display 
of  vertical- 
component  Lg 
signals  for 
the  current 
event. 


the  piupii^iiiion  plain  i,  ii'  ^  ^  .  ..j  ,  .  .:_ii  uicse 

two  points. 


The  seismic  data  themselves  can  iv  accessed  hy  relumine  i  >  the 
main  menu  (Figure  2)  and  sequentially  selecting  a  seismic  phase  from  the 
PHASE  menu  and  Analysi  Station  from  the  FUNCTION  menu.  'Fhus,  for 
example,  specifying  the  phase  1.^  in  tins  maiUicr  auioinalically  initiates  a 
process  by  which  the  L,,  recordings  liom  ihe  JVH  explosion  are  extracted 
from  the  on-line  databa.se  of  digital  waveform  data  and  displayed  on  the 
screen  in  the  format  shown  in  Figure  11.  In  this  and  subsequent  seismo¬ 
gram  displays,  the  stations  arc  ordered  bv  increasing  epicentral  distance 
and,  for  the  Lg  and  .surface  waves,  the  data  aie  plotted  as  a  function  of 
group  velocity.  It  can  be  .seen  in  this  figure  that  the  traces  are  marked  by  a 
vertical  line  at  a  group  velocity  of  3,5  km/sec,  indicating  the  nominal 
expected  onset  time  of  the  Lg  phase.  This  display  permits  the  analyst  to 
quickly  assess  the  quality  of  the  data  and  identify  any  prominent  character¬ 
istics.  Thus,  for  example,  it  can  ire  .seen  from  Figure  1 1  that  the  Lg  signals 
are  relatively  weak  at  the  IRIS  station  KIV  (IKIV)  and  at  NORSAR 
(NOIAO),  consistent  with  liic  expe.. teii  efh  k  (>f  propagation  across  the 
Ca.spian  Sea  (KIV)  anil  along  the  complex,  fa  legional  path  to  Norway. 
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RMSLg 

jtartup 


Figure  12. 

liagnitude 

Measurement 

nenu. 


Tools 


Having  vcriliod  that  the  available  Lg  data  are  suitable  tor  further 
pi  Ke-  'he  analyst  can  proceed  to  magnitude  estimation  by  selecting 
the  Ma^’nitiide  Measurement  option  from  the  FUNCTION  menu  as  illus- . 
trated  in  Figure  12.  Selecting  RMS  Lg  from  this  li.sl  of  live  currently  avail¬ 
able  magnitude  measures  initiates  a  scries  of  operations  by  which  the  Lg 
wa'  eforms  recorded  from  the  JVE  event  are  automatically  processed  using 
algorithms  described  by  Ringdal  (1983)  to  obtain  single-station  and  net- 
work-averaged  measures  of  the  Lg  magnitude.  Mj  g.  Since  this  processing 
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lakes  some  time  to  accompli.sh,  it  is  run  in  the  background,  thereby  permit¬ 
ting  the  analyst  to  move  on  to  other  tasks  while  it  is  being  completed.  In  the 
present  example,  the  analy.st  has  selected  Teleseismic  P  from  the  PHASE 
menu  and  Analyst  Station  from  the  FUNCTION  menu,  resulting  in  the 
waveform  display  shown  in  Figure  13.  The  P  wave  seismograms  are  plot¬ 
ted  here  as  a  function  of  reduced  lime,  extending  from  10  seconds  before  to 
15  seconds  after  the  signal  on.set  time  which  is  denoted  by  the  vertical  line 
segmenus.  Generally,  there  are  more  recordings  than  can  be  displayed  on  a 
single  .screen  in  this  formal  and,  in  such  ca.scs,  the  analyst  can  readily  page 
forward  and  backward  through  these  multiple  screens  using  the  NEXT  and 
PREVIOUS  menu  buttons  at  the  top  of  the  analy.st  station  display. 


An  extensive  .set  of  signal  priKCs.sing  and  analysis  capabilities  has 
been  assembled  within  the  analy.st  station  module,  as  indicated  by  the 
TOOLS  menu  display  in  Figure  1 3.  The.se  tools  permit  the  analyst  to; 


redefine  anivai  times  {Pick  Arrival ), 
a.s.sign  waveform  quality  [Set  Quality), 
measure  amplitudes  and  periods  (Expand,  m/,). 
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Filter 


Compare 


Figure  13. 

Analyst  sta¬ 
tion  display 
of  selected 
vertical-com¬ 
ponent 
teleseismic  P 
wave  signals 
for  the  current 
event. 


•  filter  the  daia  (Demf.yyuvi  j,. 

•  normalize  to  a  common  instrument  respc.ise  t /rs."  >') 

•  compare  a  selected  wavetorm  will,  wii'/etoini.'  ict-c!  le  1  n  .  .-j;  a- 

tion  from  previous  explosions  (Compare, Clone). 

Thus,  for  example,  it  can  be  seen  from  Figure  1 3  that  the  P  wave 
signal  recorded  at  station  MAJO  in  Japan  (bottom  a  scured  by 

microseismic  noise  lying  outside  the  signal  passband.  Selecting  filter  from 
the  TOOLS  menu  causes  the  interactive  filler  interface  to  be  displayed  as 
shown  in  Figure  14,  where  the  analyst  has  used  the  mouse  to  adjust  the 
slider  bars  to  define  low  and  high  frequency  cutoffs  of  0.5  and  5.0  Hz, 
respectively,  and  to  specify  a  second  order  roll-off  outside  that  band  by  set¬ 
ting  the  corresponding  pole  numbers  to  2.  Subsequent  selection  of  any 
trace  with  the  mouse  causes  this  filter  to  be  automatically  applied  and  the 
result  displayed,  as  shown  here  for  station  MAJO. 

Application  of  the  Compare  function  from  the  TOOLS  menu  pro¬ 
vides  another  illustration  of  the  powerful  capability  foi  interactive  analysis 
which  is  available  in  the  analyst  station.  For  example,  by  sequentially 
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Figure  14. 

.example  of 
the  specifi¬ 
cation  and 
subsequent 
application 
of  a  band¬ 
pass  filter  to 
the  data  of 
Figure  13. 


selecting  this  function  and  the  station  KONO  waveform  from  Figure  14,  a 
search  through  the  database  is  automatically  initiated  to  identify  all  other 
available  P  waveforms  recorded  at  that  station  from  previous  Shagan  River 
explosions.  Generally,  not  all  of  the  explosions  identified  by  this  search  are 
of  equal  importance  for  comparative  purpo.scs  and  selection  of  a  meaning¬ 
ful  subset  is  facilitated  by  examining  their  locations  in  the  context  of  the 
SPOT  image  te.st  site  information  interface,  as  indicated  in  Figure  15. 
Using  this  interface,  the  analy.st  can  interactively  .select  specific  events  with 
the  mouse,  as  indicated  here  by  the  diamond  overlays,  on  the  basis  of  prox¬ 
imity  to  the  current  event  or  .some  other  criterion,  and  display  the  corre¬ 
sponding  P  waveforms  in  a  time  expanded  analyst  station  mode  such  as 
that  of  Figure  16.  Here  the  current  JVE  event  recording  is  shown  at  the  top 
of  the  figure  and  the  waveforms  from  the  selected  comparison  events  are 
displayed  beneath  it  in  the  order  in  which  they  were  .selected. 

Cions  Also  illustrated  in  Figure  16  is  the  u.se  of  the  Clone  feature  from  the 

TOt  )LS  menu  which  permits  the  operator  to  .select  any  trace  on  the  screen 
with  the  mouse  (in  this  case  the  JVE  recording),  create  a  color-coded  (red) 
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Figure  15. 

SPOT  locations 
of  Shagan  River 
explosions 
recorded  at  sta¬ 
tion  KONO  in 
Norway.  The  dia¬ 
mond  symbols 
denote  those 
events  selected 
for  comparative 
analysis. 


copy  and  drag  it  to  any  other  location  on  the  screen  to  conduct  detailed 
waveform  comparisons.  This  capability  permits  the  analyst  to  quickly 
assess  the  consistency  of  the  current  observation  with  past  experience  and, 
in  this  case,  to  conclude  that  the  JVE  recording  at  KONO  is  quite  consis¬ 
tent  with  previous  oh.scrvaiions  at  that  station  iVoin  other  explosions 
located  in  that  area  of  the  test  site. 


Figure  16. 

Comparison  of 
the  P  wave  sig¬ 
nal  recorded  at 
the  station 
KONO  from  the 
current  event 
(top  and  red) 
with  the  signals 
recorded  at  that 
station  from  the 
selected  events 
of  Figure  15. 


Expand 
and  mi, 


Figure  17. 

Illustration  of 
the  interactive 
determination 
of  the  ampli¬ 
tude  and 
period  values 
used  to  define 
single  station 
mf,  values. 
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o.i.'  hiitiai  waveiorm  r-\  u  tias  tin'ii  completed,  the  analyst 
c;  !i  p!>  .  J  i'>  majtnitnde  estimation  nsine  llie  Expand  and  n./,  (unctions 
trom  the  IDOLS  menu.  By  sequentially  selecting  Expand  and  any  wave¬ 
form  from  the  analyst  station  display,  a  magnilied  reproduction  of  that 
waveform  is  produced  as  shown  in  Figure  17.  Then,  by  employing  the 
lunction,  the  mou.se  can  be  used  to  position  a  rectangle  so  as  to  define  the 
peak-to-peak  amplitude  and  half  period  of  the  selected  cycle  of  motion  as 
shown  in  this  figure.  After  this  proce.ss  has  been  repeated  for  each  usable 
trace,  the  resulting  amplitudes  are  input  to  the  nij,  estimation  module  using 
the  Magnitude  Measurement  menu  of  Figure  12  where  they  are  converted 
to  ground  motion,  corrected  for  epicentral  distance  and  station  effects  and 
logarithmically  averaged  to  obtain  a  network-averaged  m(,  value  and  asso¬ 
ciated  uncertainty.  The  resulting  individual  station  m^,  values  for  the  cur¬ 
rent  JVE  explosion  are  displayed  in  Figure  18,  together  with  the  estimated 
network-averaged  value  of  6.012.  This  display  module  also  provides  the 
capability  for  the  analyst  to  interactively  eliminate  questionable  data  points 
and  then  to  recalculate  the  network  average  using  the  ADD/DELETE  and 
RECALC  menu  buttons  shown  at  the  to*'  ~ ''  the  figure.  This  process  can  be 
continued  until  a  final  stable  e^^imate  is  obtained,  at  which  time  it  is 
written  to  the  database  using  the  designated  SAVE  menu  button. 
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Figure  18. 

Comparison  of 
individual  sta¬ 
tion  and  net¬ 
work-averaged 
nif,  magnitudes 
determined  for 
the  current 
event. 


RMS  Lg  Note  from  Figure  18  that  at  this  point  the  system  has  notified  the 

analyst  that  the  RMS  Lg  magnitude  estimation  has  been  completed,  as  indi¬ 
cated  by  the  icon  in  the  upper  left  hand  coinei  of  tins  display.  Selecting  this 
icon  with  the  mouse  produces  the  display  of  individual  station  and  net¬ 
work-averaged  M|  g  values  shown  in  Figure  19.  It  can  t.e  .seen  from  this  fig¬ 
ure  that  the  estimated  network-averaged  Mj  g  value  of  5.967  is  somewhat 
lower  than  the  coiresponding  mj,  value  of  6.012  of  Figure  18,  and  it  would 
be  natural  for  the  analyst  to  question  whether  there  is  any  significance  to 


Figure  19. 

Comparison  of 
individual  sta¬ 
tion  and  network- 
averaged  Mig 
magnitudes 
determined  for 
the  current  event. 


IS 


Figure  20. 

Comparison  of 
the  observed 
value  ofmf,- 
(NOR- 
SAR)  for  the 
current  event 
with  the  con¬ 
tours  represent¬ 
ing  observed 
variation  of 
that  parameter 
for  previous 
events. 

New 

Instrument 
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tins  JiiK'ience  This  issue  can  he  addressed  by  selecting  a  representative 
stutir:  V  g.  N't)Rr.)  front  the  RU.SIDUAL  menu  on  this  display  and  return¬ 
ing  to  the  SPOT  image  lest  site  information  interface  to  overlay  the  quan¬ 
tity  mj,  -  Mj  g(N()RS)  at  the  current  event  location  as  in  Figure  20,  where 
contours  of  this  quantity  derived  from  results  of  previous  explosions  have 
also  been  overlaid  for  reference  purposes.  It  can  be  seen  from  this  display 
that  while  this  magnitude  difference  is  typically  negative  in  the  northeast 
quadrant  of  the  lest  site  (by  as  much  as  0.15  magnitude  units),  it  is  gener¬ 
ally  positive  in  the  vicinity  of  the  current  JVE  event  location,  with  an  aver¬ 
age  value  close  to  the  observed  value  of  0.043.  Thus,  this  capability 
permits  the  analyst  to  quickly  conclude  that  the  current  result  is  entirely 
consistent  with  previous  experience  in  this  area  of  the  test  site. 
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The  long-period  surface  wave  data  can  be  accessed  in  the  same 
manner  as  the  Lj,  and  P  pha.se  data  described  above,  resulting  in  analyst  sta¬ 
tion  displays  such  as  that  shown  in  Figure  21  for  the  vertical  component 
Rayleigh  wave  phase.  It  can  be  seen  from  this  example  that,  although  these 
o'ices  are  aligned  according  to  a  common  group  velocity  scale,  the  surface 
wave  signals  are  quite  viu-iable  and  difficult  to  correlate  from  trace  to  trace. 
Tl  ’-  is  at  least  partially  due  to  the  fact  that  the  re.sponse  characteristics  of 
the  long-period  instruments  used  to  record  these  data  vary  significantly 
between  stations.  Thus,  for  example,  the  station  MAJO  trace  shown  at  the 
bottom  of  this  figure  was  recorded  through  a  relatively  broadband  system. 
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Figure  21. 

Analyst  sta¬ 
tion  display 
of  selected 
long-period 
Rayleigh 
wave  signals 
for  the  current 
event. 


Figure  22. 

Example  of  the 
specification 
and  subsequent 
application  of 
the  instrument 
response  nor¬ 
malization  fea¬ 
ture  to  the  data 
of  Figure  21. 


with  the  result  that  the  long-period  signals  of  principal  interest  are 
obscured  by  higher  frequency  arrivals.  In  order  to  provide  a  common  basis 
for  comparison,  the  TOOLS  menu  includes  a  function  (New  instrument) 
which  permits  the  analyst  to  transform  the  data  to  that  which  would  have 
been  observed  if  the  same  instrumentation  had  been  employed  at  each  sta¬ 
tion.  This  feature  is  graphically  illusU-ated  in  Figure  22  where  the  analyst 
has  selected  the  nominal  SRO  long-period  response  and  converted  all  of 
the  naces  in  this  display  to  that  response  by  simply  designating  them  with 
the  mouse.  It  is  evident  from  this  example  that  such  a  capability  greatly 
facilitates  quantitative  comparisons  and  evaluation  of  data. 
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SpBCtTBl  P  Once  ific  available  lung-pcnou  Rayicigh  and  Love  wave  data  have 

boon  r'  oviewcd,  the  analyst  can  return  to  the  Magnitude  Measurement 
menu  ol  Figure  (2  and  select  Moment  Tensor,  (hereby  initiating  a  process 
by  which  the  observed  long-period  data  are  phase  matched  filtered  and  for¬ 
mally  inverted  to  obtain  a  surface  wave  measure  of  seismic  magnitude. 
While  this  computationally  intensive  process  is  running  in  the  background, 
the  analyst  next  selects  Spectral  P  from  this  same  menu  and  proceeds  with 
an  interactive  determination  of  the  network-averaged  P  wave  spectrum  cor¬ 
responding  to  the  current  event.  The  spectrum  estimated  in  this  process  can 
be  inverted  to  obtain  either  an  equivalent  spectral  magnitude  or  a  direct, 
model-based  yield  estimate  (Murphy,  1989),  and  Figure  23  shows  the  menu 
structure  for  the  latter  option,  where  it  is  indicated  that  the  Mueller/Murphy 
granite  (M/M  Granite)  source  model  has  been  selected  for  the  current  JVE 
event.  The  resulting  automatic  fit  to  the  attenuation-corrected,  observed 
spectrum  is  shown  in  Figure  24  where  it  is  indicated  that  a  theoretical  spec¬ 
trum  corresponding  to  a  yield  (W)  of  119  kt,  a  pP  -  P  delay  time  (T^)  of 
0.82  seconds  and  a  pP/P  amplitude  ratio  (A)  of  0.17  provides  the  best  over¬ 
all  fit  to  the  data.  Alternately,  the  analyst  can  elect  to  interactively  change 
these  automatically  determined  model  parameters  using  the  plus  and  minus 
buttons  located  in  the  box  in  the  upper  right  corner  of  this  display,  to 
explore  solutions  corresponding  to  goodness  of  fit  criteria  not  incorporated 
in  the  automatic  algorithm. 

Figure  23. 

Menu  options 
for  the  estima¬ 
tion  of  network- 
averaged  P 
wave  spectra. 
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Figure  24. 

Comparison  of 
normalized 
observed  and 
best-fitting  theo¬ 
retical  network- 
averaged  P 
wave  spectra 
for  the  current 
event. 


The  icon  appearing  in  the  upper  left  corner  of  Figure  24  indicates 
that  the  surface  wave  phase  matcheJ  filtering  processing  has  now  been 
completed,  and  selecting  it  with  the  mouse  provides  access  to  graphical 
displays  of  these  processing  results.  Following  review  of  these  results,  the 
operator  can  select  the  associated  MTIS  menu  option  to  view  a  graphical 
representation  of  the  resulting  moment  tensor  solution  shown  in  Figure  25 
where  the  model  fits  to  the  Rayleigh  (left)  and  Love  (right)  wave  data  cor¬ 
responding  to  a  tectonic  release  F  factor  of  0.13  are  displayed.  It  is  infor¬ 
mative  to  examine  this  .solution  in  the  context  of  other  available 
information  regarding  the  te.st  environment.  This  can  be  accomplished  by 


Comparison  of  the 
inferred  surface 
wave  moment  ten¬ 
sor  solution  for 
current  event  with 
the  corresponding 
path  normalized, 
observed  Rayleigh 
(left)  and  Love 
( right)  wave 
amplitude  data. 


Moment 

Tensor 


Figure  25. 
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Figure  26. 

Comparison  of 
the  surface  wave 
moment  tensor 
solution  for  the 
current  event 
(yellow  and  light 
blue  concentric 
circles)  with 
those  for  nearby 
Shagan  River 
explosions  (red 
and  dark  blue 
concentric  cir¬ 
cles)  and  with  the 
surface  geologic 
map  of  the  area. 


lmajt:-^ijLccatlonsi 


CK'erlays 


IRefresW 


iF unction 


Help 


iQuitj 


returning  once  again  to  the  SPOT  image  test  site  information  interface  and 
successively  overlaying  the  current  and  historical  moment  tensor  solutions, 
as  welt  as  the  map  of  surface  geologic  features,  as  shown  in  Figure  26.  In 
this  display,  the  concentric  circles  represent  the  moment  tensor  solutions 
where,  for  the  historical  events,  the  ratio  of  the  diameters  of  the  red-to-blue 
circles  is  equal  to  the  inferred  tectonic  release  F  factors  and  the  line  seg¬ 
ments  interior  to  the  circles  are  parallel  to  the  associated  strike  of  the  tec¬ 
tonic  component.  The  solution  for  the  current  (JVE)  event  is  differentiated 
here  by  an  alternate  color  .scheme  in  which  yellow  and  light  blue  circles 
replace  red  and  blue,  respectively.  It  can  be  seen  from  this  display  that  the 
.solution  for  the  current  event  is  quite  con.sistent  with  those  previously 
determined  for  nearby  historical  events.  Moreover,  it  is  evident  that  the 
inferred  strikes  of  the  tectonic  relea,scs  triggered  by  these  explosions  are 
parallel  or  subparallel  to  that  of  the  mapped  Chinrau  fault.  Thus,  by  using 
the  test  site  information  interface,  the  analyst  is  able  to  quickly  verify  that 
the  surface  wave  moment  ten.sor  solution  which  was  obtained  for  the  cur¬ 
rent  event  is  consistent  with  both  the  .solutions  from  previous  nearby  explo¬ 
sions  and  with  the  regional  tectonic  environment  in  which  the  test  was 
conducted. 
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Yield  Esth 
mation 


Having  determined  the  various  seismic  measures  of  source  size 
(including  a  regional  magnitude  not  discussed  above),  the  analyst  can 
proceed  to  estimate  the  corresponding  explosion  yield  and  the  associated 
uncertainty  in  that  estimate  by  selecting  the  Yield  Estimation  feature  from 
the  FUNCTION  menu  of  Figure  2.  This  yield  estimation  module  provides 
an  estimate  of  unified  yield  based  on  multiple  magnitude  measures  and 
associated  extremal  confidence  limits  corresponding  to  different  sets  of 
constraints  specified  by  the  analyst.  The  confidence  limits  are  obtained  as 
the  solution  of  a  nonlinear  programming  problem  in  which  estimates  of 
yield  for  a  training  set  of  explosions  are  minimized  and  maximized  over 
the  space  of  admissible  parameters  (Rodi,  1989;  Rodi  and  Murphy,  1990). 
Input  to  this  model  consists  of  multiple  network-averaged  magnitudes  for 
the  current  event  and  a  nominal  m^/yield  relation  based  on  some  combina¬ 
tion  of  data  analysis  and  expert  opinion,  which  for  explosions  at  the  Shagan 
River  test  site  is  taken  to  be  (Murphy,  1990); 


=  4.45  +■  0.75  logW 


Figure  27. 

Menu  and 
sample  out¬ 
put  for  the 
unified  yield 
estimation 
module. 


The  current  module  can  provide  estimates  corresponding  to  five 
different  sets  of  constraints,  parameterized  by  the  uncertainty  in  the  aver¬ 
age  value  of  mi5  at  100  kt  (m|,(100)),  the  uncertainly  in  the  slope  (b)  of  the 
mj/yield  relation  and  the  upper  bound  on  the  absolute  value  of  the  correla- 
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Flgue  28. 

Comparison  of 
unified  yield 
estimates  (W) 
and  associated 
uncertainties 
( F)  obtained 
with  (top)  and 
without  (bot¬ 
tom)  the  sur¬ 
face  wave 
moment  tensor 
magnitude. 


lion  oilicicnu  uinongst  ihe  yield  cViiiiiaiion  errors  lor  Uie  individual  mag- 
niuidos  'max  P).  The  menu  for  the  yield  estimation  module  is  shown  in 
Figure  27  together  with  sample  output  corresponding  to  the  application  of 
Models  1  and  3  to  the  magnitude  data  for  the  current  JVE  event.  In  these 
displays,  the  resulting  unified  yield  estimates  (W)  and  associated  upper 
bound  uncertainty  factors  (F)  are  .shown,  together  with  the  yield  estimates 
obtained  from  the  individual  magnitude  values.  Note  that  the  yield  value 
is  shown  here  as  an  open  circle  to  denote  the  fact  that  it  was  not  used  in  the 
formal  computation  of  W  and  F  due  to  the  fact  that  an  adequate  training  set 
of  previous  values  does  not  currently  exist  for  that  magnitude  measure.  It 
can  be  seen  that  for  these  two  ca.ses,  the  unified  yield  estimate  is  about  121 
kt  and  the  uncertainty  bounds  range  from  a  factor  of  1.96  for  the  case 
(Model  1)  in  which  mj,{100)  is  specified  as  5.95  ±  0.05,  to  a  factor  of  1.66 
for  the  case  (Model  3)  in  which  mj„(  100)  is  given  as  exactly  5.95.  In  Figure 
28,  the  Model  1  .solution  from  Figure  27  (top)  is  compared  with  the  result 
of  rerunning  Model  1  without  the  surface  wave  moment  tensor  magnitude 
(i.e.,  by  using  the  slider  bars  at  the  bottom  of  the  menu  box  to  assign  a 
weight  of  zero  to  W(Mo)).  It  can  be  seen  that  this  modification  results  in  a 
significant  reduction  of  the  F  factor  from  1.96  to  1.55,  reflecting  the  rela¬ 
tively  uncertain  yield  estimation  capability  of  the  surface  wave  moment  as 
determined  from  the  analysis  of  the  training  set  of  multiple  magnitude  data. 
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Statistical 

Summary 


Given  an  optimum  seismic  estimate  of  unified  yield  and  the  associ¬ 
ated  uncertainty,  the  analyst  can  next  proceed  to  assess  whether  this  esti¬ 
mate  is  consistent  with  the  150  kt  threshold  of  the  Threshold  Test  Ban 
Treaty  (TTBT)  by  selecting  the  Statistical  Summary  option  from  the 
FUNCTION  menu.  This  module  provides  access  to  a  variety  of  statistical 
compliance  assessment  tests  developed  by  Mission  Research  Corporation 
(Gray  etal.,  1990),  as  indicated  by  the  menu  displayed  in  Figure  2  ’.  In  this 
example,  the  analyst  has  requested  a  summary  of  the  results  of  three  single 


Figure  29. 

Statistical 

assessment 

menu. 


event  compliance  tests,  producing  the  graphical  summary  shown  in  Figure 
30,  where  the  red  areas  under  the  distribution  curves  provide  a  measure  of 
the  probability  of  a  violation  of  the  150  ki  threshold.  The  formal  statistical 
results  are  summarized  in  the  boxes  in  the  upper  left  hand  corners  of  these 
figures,  and  it  can  be  seen  that  the  results  of  each  of  these  three  tests  indi¬ 
cate  that  the  seismic  data  arc  consistent  with  a  yield  of  less  than  150  kt  at  a 
specified  false  alarm  rate  of  less  than  2.5  percent.  The  resolving  power  of 
such  tests  can  be  illustrated  by  considering  the  outcome  which  would  have 
resulted  if  these  same  seismic  data  had  been  observed  from  a  test  below  the 
water  table  at  NTS.  In  this  case,  the  Shagan  mj/yield  relation  would  be 
replaced  by  (Murphy,  1981); 

nijj  =  3. +  0.81  lof’W 


26 


(hen  iew  oj  Sele  led  YFS  Capabilities 


Figure  30. 

Comparison 
of  the  results 
of  three  differ¬ 
ent  tests  of 
seismic  com¬ 
pliance  of  the 
current  event 
with  the  150 
kt  threshold  of 
the  TTBT 
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resulting  in  the  revised  Test  1  outcome  shown  in  Figure  31.  It  can  be  seen 
that  under  this  hypothesis,  the  yield  estimate  would  be  362  kt  and  it  would 
be  concluded  that  this  explosion  violated  the  150  kt  limit  of  the  TTBT. 
Thus,  this  simple  example  graphically  illustrates  the  importance  of  the  test 
site  magnitude  bias  effect  which  has  received  such  intense  study  since  the 
initial  negotiation  of  the  TTBT. 


Figure  31. 

Compliance 
test  results  (Test 
I)  for  the  sce¬ 
nario  in  which 
the  current 
event  data 
were  observed 
from  an  explo¬ 
sion  below  the 
water  table  at 
NTS. 


27 


Overview  of  Selected  YES  Capabilities 


Spread¬ 

sheet 


Having  completed  the  formal  analysis,  the  Spreadsheet  option  from  the 
FUNCTION  menu  can  be  used  to  provide  a  less  formal  environment  for  evalu¬ 
ating  the  seismic  yield  estimate  for  the  current  event  in  the  context  of  the 
results  obtained  from  analyses  of  previous  explosions  at  that  test  site,  as  well  as 
any  available  independent  yield  calibration  data,  as  illustrated  in  Figure  32.  In 
this  display,  the  four  individual  seismic  yield  estimates  and  corresponding  uni¬ 
fied  yields  are  listed  for  forty  of  the  largest  Shagan  River  explosions,  together 
with  selected  calibration  yields  (W (Bocharov)),  taken  here  for  the  purposes  of 
illustration  to  be  the  three  high  yield  values  published  by  Bocharov  et  al. 
(1988).  The  ratios  of  the  unified  seismic  to  calibration  yields  are  listed  in  the 
last  column,  together  with  the  average  offset  (Avg  Bias)  and  uncertainty  (F) 
inferred  from  these  limited  calibration  data.  For  the  nominal  magnitude/yield 
relations  listed  at  the  top  of  this  display,  the  average  offset  is  a  factor  of  1.07 
with  an  associated  F  factor  of  1.25.  Figure  33  shows  the  result  of  interactively 
adjusting  each  magnitude/yield  relation  in  such  a  manner  that  the  individual 
average  offsets  are  minimized.  It  can  be  seen  that  these  modifications  eliminate 
the  average  bias  in  the  unified  seismic  yield  estimate  (i.e.,  Avg  Bias  =  1.00) 
and  reduces  the  F  factor  to  1.20.  This  example  illustrates  how  the  spreadsheet 
module  provides  a  capability  for  the  analyst  to  rapidly  explore  the  effects  of 
alternate  hypotheses  on  the  seismic  yield  estimation  process. 


Figure  32. 

Spreadsheet 
summary 
comparison 
of  seismic 
yield  esti¬ 
mates  for  the 
current  event 
(09/14/88) 
with  those 
obtained  for 
selected  pre¬ 
vious  Shagan 
River  explo¬ 
sions. 
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Figure  33. 

Spreadsheet 
summary 
illustrating 
the  results  of 
interactively- 
modifying  the 
designated 
magnitude/ 
yield  rela¬ 
tions. 


Event 

Summary 
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Once  the  procejssing  has  been  completed  for  the  current  event,  the 
analyst  can  produce  a  one-page  summary  of  the  analysis  results  by  select¬ 
ing  the  Event  Summary  option  from  the  FUNCTION  menu.  As  is  indicated 
in  Figure  34,  listed  results  include  the  date,  origin  time  and  location  of  the 
explosion,  the  derived  unified  seismic  yield  and  its  associated  uncertainty, 
the  results  of  the  statistical  test  of  compliance  of  the  current  event  with  the 
150  kt  limit  of  the  TTBT,  a  tectonic  release  characterization  of  the  event 
based  on  the  surface  wave  moment  tensor  analysis  and,  if  available,  the 
CORRTEX  yield  estimate  for  the  explosion.  In  addition,  each  of  the  esti¬ 
mated  seismic  magnitudes  is  listed  together  with  its  associated  yield  esti¬ 
mate  and  assigned  weight  in  the  determination  of  the  unified  yield  estimate. 


The  sample  analysis  .sc.ssion  described  above  has  provided  an  over¬ 
view  of  some  of  the  capabilities  which  are  currently  available  within  the 
Yield  Estimation  System.  As  was  evident  from  the  menu  structures  shown 
in  the  graphical  displays,  there  are,  in  addition,  many  other  features  which 
were  not  exercised  here  in  order  to  hold  the  description  to  a  manageable 
length.  A  complete  description  of  all  currently  available  system  options  is 
provided  in  Appendix  C. 
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Overview’  of  Selected  YES  Capabilities 


Figure  34. 

Summary  of 
analysis  results 
for  the  selected 
event. 
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Summary 


In  this  report  we  have  presented  an  in-depth  description  of  the  cur¬ 
rent  status  of  the  on-going  research  investigations  directed  toward  the 
development  of  a  comprehensive  new  seismic  yield  estimation  system  for 
underground  nuclear  explosions.  In  its  current  configuration,  the  YES  sys¬ 
tem  is  applicable  to  underground  nuclear  explosions  at  the  Soviet  Shagan 
River  and  Novaya  Zemlya  test  sites  and  encompasses  a  database  of  more 
than  10,000  digital  seismograms  recorded  at  stations  of  the  USAEDS, 
GDSN,  CDSN  and  IRIS  networks  from  explosions  at  these  two  sites.  For 
both  test  areas,  information  regarding  the  explosion  source  environment  is 
presented  to  the  analyst  in  the  context  of  SPOT  satellite  images  of  the  sites, 
together  with  associated  surface  and  subsurface  geologic  information  and 
DMA  topographic  data.  The  on-line  database  for  YES  also  contains  a  wide 
variety  of  tabular  information,  including  complete  event  and  station  loca¬ 
tion  files  containing  both  classified  and  unclassified  locations,  standard 
travel-time  tables  for  the  seismic  arrivals  used  for  yield  estimation,  propa¬ 
gation  path  and  station  corrections  for  use  in  magnitude  determinations  and 
a  comprehensive  instrument  response  database.  The  instrument  database 
encompasses  a  compilation  of  over  750  different  response  histories  col¬ 
lected  from  the  USGS,  AFTAC,  IRIS  and  others  which  detail  the  character¬ 
istics  of  the  various  recording  instruments  as  a  function  of  station,  channel 
and  date. 

The  capabilities  and  functionality  of  the  current  version  of  the  YES 
system  were  graphically  illustrated  in  Section  11  using  displays  of  the 
.screens  encountered  by  an  analyst  in  a  typical  processing  session.  The 
selected  capabilities  displayed  in  this  sample  session  were  related  to  a  more 
complete  definition  of  the  system  functionality  through  references  to  a 
series  of  appendices  containing  a  detailed  script  of  the  operator  actions 
required  to  reproduce  the  sample  session  (Appendix  B),  a  complete 
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Summary 


description  of  all  currently  available  system  options  (Appendix  C),  as  well  as 
information  regarding  software  and  hardware  implementation  requirements 
(Appendix  A)  and  a  top  level  schematic  diagram  of  the  main  modules  of  the 
system  and  their  interdependencies  (Appendix  D).  Taken  together.  Section  IT 
and  its  associated  appendices  provide  a  concise  description  of  the  current  sta¬ 
tus  of  the  YES  system  development  effort. 
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Requirements  and  Startup 
Instructions 


Hardware  Requirements 

The  YES  system  is  designed  to  run  on  SUN-4  and  SUN  SPARCSta- 
tion  computers  with  a  minimum  of  32MB  of  resident  memory.  It  uses  the 
standard  SUN  8-bit  color  monitor  with  a  resolution  of  1152x900  pixels,  as 
well  as  the  standard  SUN  3-button  mouse.  YES  also  requires  a  minimum  of 
I  gbyte  of  disk  storage  to  hold  the  complete  system  and  its  required  data. 

Software  Requirements 

YES  is  written  in  C  and  FORTRAN  and  uses  the  X  Window  Sys¬ 
tem  software,  version  X 1 1  revision  4.  It  also  uses  the  Hewlett  Packard  wid¬ 
get  set  extension  to  X  and  the  Motif  widget  set,  release  1.1.  It  has  been 
successfully  compiled  using  the  standard  SUN  bundled  C  compiler  and  the 
SUN  1.3.1  FORTRAN  compiler,  and  is  designed  to  run  under  the  SUN 
operating  system  version  SUN-OS  4. 1 .0  or  later. 

Startup  instructions 

The  following  instructions  assume  that  the  X  Window  System  and 
YES  have  been  installed  according  to  S-CUBED  specifications.  Also  refer 
to  the  information  on  mouse  and  window  techniques  in  Appendix  C  before 
proceeding. 

To  bring  up  the  X  environment  and  run  YES  you  must  do  the  fol¬ 
lowing: 

I .  Login  to  UNIX  at  the  console. 
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2.  Add  /usr/bin/Xll  to  your  search  path,  if  it  is  not  already  there,  and 
issue  the  command: 

prompt>  rehash 

(where  prompt>  denotes  the  UNIX  operating  system  prompt). 

3.  Ensure  that  the  hies  xinitrc  and  twmrc  as  provided  on  the  distribu¬ 
tion  tape  have  been  copied  to  .xinitrc  and  .twmrc  in  your  home 
directory. 

4.  Issue  the  command: 

prompt>  startx 

and  wait  until  an  X  terminal  window  appears  on  your  screen. 

5.  Move  the  cursor  into  the  X  terminal  window  and  change  to  the 
main  directory  where  YES  has  been  installed. 

6.  Issue  the  command: 

proinpt>  source  Init 

7.  Issue  the  command: 

proinpt>  yieldsys& 


Shutdown  Instructions 

To  shutdown  the  system: 

1.  Choose  Quit  from  the  YES  main  menu. 

2.  Exit  X  by  pointing  outside  any  window.  Click  and  hold  the  right 
button.  Move  the  pointer  down  to  menu  item  ExitX  Windows,  stop, 
then  move  to  the  right  and  release  on  Select  if  you  're  sure. 

3.  Logout  of  UNIX  at  the  con.snle. 
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Script  For  Sample  Session 


In  this  section,  all  of  the  operator  actions  required  to  duplicate  the 
sample  session  shown  in  Section  II  are  specified.  Note  that  when  run  in  a 
classified  environment,  the  values  of  some  of  the  parameters  and  process¬ 
ing  results  will  be  somewhat  different  from  those  shown  here.  A  full 
description  of  all  system  options  is  provided  in  Appendix  C,  including 
details  regarding  menu  structures  and  the  mouse  actions  required  to  initiate 
specific  actions.  All  figures  mentioned  in  this  section  are  found  in  Chapter 
2  of  this  manual. 

The  procedures  required  to  initiate  and  start  the  program  are  speci¬ 
fied  in  Appendix  A.  To  recap: 

1.  Change  to  the  YES  directory. 

2.  T>rpe  source  Init. 

3.  yieldsys&. 

Once  the  program  has  been  started,  the  first  screen  display  is  the 
title  page.  Clicking  anywhere  within  the  display  brings  up  the  main  menu. 
A  site  and  an  event  are  normally  selected  as  the  first  step  in  the  analysis. 

1.  Choose  Balapan  (Shagan)  from  the  SITE  menu. 

2.  Choose  Shagan  Events->I980-I989->I988->September  14  (JVE) 
from  the  EVENT  menu. 

This  event  will  be  referred  to  as  the  current  event. 

After  a  site  and  event  have  been  chosen,  one  can  proceed  to  analyze 
the  near-source  environment  of  the  current  event. 


Mam 

Menu 


Starting 
the  pro¬ 
gram 
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1.  Choose  SateUUe  unage  from  the  FUNCTION  menu. 

This  brings  up  the  satellite  image  of  the  test  site.  It  will 
come  up  with  symbols  showing  the  locations  of  all  the  events  that  occured 
at  the  test  site  and  the  current  event  location  highlighted  (see  Figure  3). 

2.  Click  on  the  current  event  diamond. 

A  description  of  the  event  is  displayed  in  the  information 
bar  below  the  menu.  Experiment  likewise  with  other  events  (i.e.  click  on 
them). 


1.  Clide  on  Full  res  in  the  side  menu  and  then  in  the  center  of  die 
crater  located  at  the  northern  end  of  the  lake. 

Once  the  10  m  resolution  image  of  the  crater  has  come  up 
(see  Figure  4).  full  resolution  images  of  adjacent  areas  can  be  viewed  by 
clicking  near  the  edges  of  the  image. 

2.  Choose  Satellite  image  from  the  IMAGES  menu. 

The  original  satellite  image  is  now  displayed. 

3.  Click  again  on  the  Full  res  button  to  turn  it  off. 

1.  Adjust  the  C  and  B  sliders  in  the  side  menu. 

This  reverses  contrast  and  alters  the  brightness  of  the  image. 

(see  Figure  5) 

2.  Click  on  the  C  and  B  buttons  to  restore  the  original  settings. 


Use  the  black  bars  to  adjust  the  sliders 


1.  Choose  Surface  geology  from  the  OVERLAYS  menu. 

2.  Choose  Current  event  from  the  OVERLAYS  menu. 

The  display  will  look  like  Figure  6. 
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Topographic  l.  Choose  Sui/ace  from  the  IMAGES  menu. 

Map  2.  Choose  Topography  contours  from  the  OVERLAYS  menu. 

3.  Choose  Current  event  from  the  OVERLAYS  menu. 

The  display  should  now  look  like  Figure  7. 

4.  Choose  Satellite  image  from  the  IMAGES  menu. 

Cross  1.  Choose  Event  locations(SPOT)  from  the  LOCATIONS  menu. 

Sactlon  2.  Click  on  Cross  sec  in  the  side  menu. 

3.  Click  below  the  current  event  location  and  again  near  the 
northern  extremity  of  the  display  (see  insert  in  Figure  8). 

The  cross  section  of  the  site  along  the  line  is  then  displayed 

(see  Figure  8). 

4.  Click  on  the  Quit  button  to  return  to  the  satellite  image. 

5.  Click  again  on  Cross  sec  in  the  side  menu  to  turn  it  off. 

Granite  l .  Choose  Depth  to  granite  from  the  IMAGES  menu. 

Surface  2.  Choose  Subsurface  contourS’>Depth  to  top  of  granite  from  the 

OVERLAYS  menu. 

3.  Choose  Current  event  from  the  OVERLAYS  menu. 

Figure  9  shows  what  the  screen  should  look  like  at  this 

point. 

4.  Choose  Satellite  image  from  the  IMAGES  menu. 

5.  Iconify  the  Satellite  Image  module. 


World 

Map 


1.  Choose  World  map  from  the  main  FUNCTION  menu. 

2.  Choose  Teleseismic  P  from  the  PHASE  menu. 

3.  Choose  Teleseismic  P->ANMO  from  the  STATION  menu. 

This  draws  the  great  circle  path  shown  in  Figure  10. 

4.  Click  on  the  REFRESH  button. 

This  removes  the  stations.  Repeat  steps  2-4  to  draw  the 
great  circle  path  for  phase  Lg  and  the  NORSAR(NOIAO)  station  shown  in 
Figure  10. 

5.  Click  on  QUIT  to  return  to  the  main  menu. 
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Analyst 

Station 


RMSLg 

Startup 


1.  Choose  Lg  from  the  PHASE  menu. 

2.  Choose  Analyst  Station  from  the  main  FUNCTION  menu. 

This  displays  the  Lg  waveforms  as  shown  in  Figure  11. 

3.  Click  on  the  QUIT  button  and  return  to  the  main  menu. 


1.  Choose  Magnitude  Measurement->RMS  Lg  from  the  FUNC¬ 
TION  menu. 

This  spawns  off  a  process,  and  when  done,  the  following 
icon  pops  up.  While  this  process  runs  in  the  background,  continue  with  the 
rest  of  the  analysis. 


The  RMS  Lg  icon 


Filter 

1.  Choose  Teleseismic  P  from  the  PHASE  menu. 

2.  Choose  Analyst  station  from  the  FUNCTION  menu. 

This  displays  the  first  six  teleseismic  P  waveforms  as  shown 

in  Figure  13. 

3.  Click  on  the  NEXT  and  PREVIOUS  buttons  to  page  through 
the  waveforms. 

4.  Choose  Filter  from  the  TOOLS  menu. 

5.  Set:  Low  Frequency  =  0.5 

Low  Pole  Number  =  2 
High  Pole  Number  =  2 


Move  the  slider  bars  or 
click  anywhere  in  the  slider 
box  to  change  the  settings 


High  Frequency  =  5.0 

6.  Click  on  the  station  MAJO  waveform. 


This  applies  the  selected  bandpass  filter  (second-order,  0.5 
to  5  Hz)  to  the  MAJO  waveform  as  shown  in  Figure  14.  To  apply  this  same 
filter  to  other  waveforms,  just  click  on  them. 

1, _ Click  on  the  EXIT  button  once  done  with  the  filter. _ 
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1.  Choose  from  the  TOOLS  menu. 

2.  Click  on  the  station  KONO  waveform. 

3.  De-iconify  the  Satellite  image  icon. 

4.  Choose  Remove  locations  from  the  REFRESH  menu  to  remove 
any  event  locations  on  the  image. 

5.  Choose  Compare  events  from  the  LOCATIONS  menu. 

6.  Click  on  Compare  in  the  side  menu. 

7.  Double-click  on  five  events  surrounding  the  current  event. 

The  display  should  look  like  Figure  1 5. 

8.  Choose  Analyst  station  from  the  Satellite  Image  FUNCTION 
menu. 

This  displays  the  KONO  Teleseismic  P  waves  for  the 
selected  events  beneath  that  for  the  JVE  event. 

1.  Choose  C/o/ie  from  the  TOOLS  menu. 

2.  Click  on  the  JVE  trace. 

A  red  copy  of  the  JVE  trace  appears  and  the  cursor  changes 
to  a  pointing  hand. 

3.  Compare  the  trace  with  the  others. 

Reposition  the  red  trace  adjacent  to  any  other  desired  U'ace 
as  illustrated  in  the  figure  above.  Also  see  Figure  16.  When  finished  with 
the  comparison,  click  again  anywhere  in  the  display  to  get  rid  of  the  clone. 


''Move  the  hand  with  the 
mouse  and  the  red  trace 
follows  it. 


4.  Click  on  die  QUIT  button  to  return  to  the  satellite  image. 

5.  Choose  Remove  locations  from  the  REFRESH  menu. 

6.  Iconify  the  Satellite  Image. 


The  Analyst  Station  display  should  now  be  on  the  screen. 
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1.  Choose  Expand  from  the  TOOLS  menu. 

2.  Click  on  station  ANTO’s  waveform. 

This  brings  up  the  waveform  in  expanded  form  as  shown  in 
Figure  17.  All  the  tools  are  still  available  for  use  in  this  mode. 

3.  Choose  ntf,  from  the  TOOLS  menu. 

4.  Measure  amplitude  and  period. 

Draw  a  rectangle  over  the  selected  half  cycle  of  motion  by 
clicking  the  cursor  at  the  peak,  dragging  the  cursor  to  the  trough  while 
holding  down  the  mouse  button  and  then  releasing  the  mouse  button  (see 
Figure  17). 

5.  Click  on  the  QUIT  button. 

The  original  full  page  screen  is  displayed.  Do  the  measure¬ 
ments  for  all  other  waveforms  using  steps  2-6.  Note  that  the  mj,  tool  can 
also  be  used  in  the  full  page  display.  Therefore,  the  measurements  need  not 
all  be  done  in  expanded  mode  although  it  is  easier  on  the  eye. 

6.  Choose  Interface  file  from  the  FUNCTION  menu. 

This  saves  the  amplitude  and  period  measurements  along 
with  other  information  for  use  by  other  modules. 

7.  Click  on  the  QUIT  button  to  get  back  to  the  main  menu. 


1.  Choose  Magnitude  Measurement~>mb  from  the  FUNCTION 
menu. 

A  display  similar  to  that  of  Figure  18  will  result.  Not  all  the 
station  magnitudes  shown  in  Figure  1 8  will  be  displayed  if  not  all  of  them 
were  measured  in  step  5  under  Expand  and  mi,. 

2.  Click  on  the  QUIT  button. 

1.  De-iconify  the  RMS  Lg  icon. 

2.  Choose  NORS  from  the  RESIDUAL  menu. 

The  display  should  now  look  like  Figure  19. 

3.  Click  on  the  QUIT  button. 

1.  De-iconify  the  Satellite  Image. 

2.  Choose  Historical  magnitude  data^>Mb-MLg(NORSAR)  from 
the  OVERLAYS  menu. 
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Instrument 

Response 


Moment  Ten¬ 
sor  Startup 


Spectral  P 


3.  Choose  Event  magnitude  data->Mb-MLg  from  the  OVERLAYS 

=’5pnii. 

The  display  should  look  like  Figure  20. 

4.  Choose  Quit  from  the  QUIT  menu. 

1 .  Choose  Rayieigh(vertical)  from  the  PHASE  menu. 

2.  Choose  Analyst  Station  from  the  FUNCTION  menu. 

The  display  should  look  like  Figure  21. 

3.  Choose  New  instrument->sro  Ip  from  the  TOOLS  menu  as 
shown  in  Figure  22. 

4.  Click  on  each  of  the  waveforms. 

This  iranstorms  all  the  waveforms  to  the  same  instrument 

re.sponse. 

5.  Click  on  the  QUIT  button. 

1.  Choose  Magnitude  Measurement->Moment  tensor  from  the 
FUNCTION  menu. 

This  spawns  a  background  process  and  the  following  icon 
pops  up  when  done.  Continue  the  analysis  while  this  is  being  processed. 

The  Xpmf  icon 


1.  Choose  Magnitude  Measurement->Spectral  P  from  the  main 

FUNCTION  menu. 

2.  Click  on  NBF. 

This  brings  up  a  display  of  the  P  wave  signal  and  noise 
spectra  for  the  first  station. 

3.  Click  on  Continue  until  all  stations  have  been  processed. 

4.  Click  on  Avgspec. 

The  resulting  display  shows  a  comparison  of  the  estimated 

network-averaged  P  wave  spectrum  (yellow  line)  with  that  expected  on 
average  for  an  explosion  having  the  current  event  m(,  value  (dashed  red 
lines). 
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Moment 

Tensor 


5.  Click  on  Yieldfit. 

6.  Click  on  M/M  Granite  (see  Figure  23) 

This  initiates  a  process  by  which  the  theoretical  solution 
corresponding  to  the  Mueller/Murphy  granite  source  model  which  provides 
the  best  fit  to  the  attenuation-corrected  observed  snert-iim  is  automatically 
determined  and  displayed.  The  yield  value  corresponding  to  this  best-fit 
solution  is  listed  at  the  top  of  the  display  together  with  its  associated  upper 
bound  uncertainty. 

7.  Click  on  Continue  twice. 

This  skips  through  the  graphs  to  obtain  the  final  aut'^  mafic 
solution  with  the  final  yield(W)  and  pP(A,To)  values  shown  in  the  box  in 
the  upper  right-hand  corner  of  the  display  (see  Figure  24). 

8.  Change  the  yield  W  from  119  kt  to  150  kt  and  then  back  down 
to  119  kt  by  using  the  +  and  -  buttons  shown  below. 


Click  on  these  +  and  -  buttons 
to  change  the  yield 


9.  Click  on  the  QUIT  button. 

This  takes  you  back  to  the  main  menu. 

1.  De-iconify  the  Xpmf  icon. 

This  brings  up  a  display  of  the  phase- matched  filter  output 
for  the  first  long-period  station.  Since  the  data  from  the  JVE  event  have 
already  been  filtered,  continue  with  the  moment  tensor  inversion  analysis. 

2.  Choose  Run  MTIS  from  the  MTIS  menu. 

3.  Click  on  some  of  the  green  and  red  squares. 

Individual  station  information  is  displayed  in  the  informa¬ 
tion  bar.  Also  see  Figure  25. 

4  Choose  SPOT  image  from  theVIEW  menu. 

This  creates  a  graphical  repre.sentation  of  the  moment  ten.sor 
solution  that  is  used  in  the  Satellite  Image  module. 
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Yield  Esti¬ 
mation 


Statistical 

Summary 


Script  for  Sample  Session 


5.  Choose  Confirm  Quit?  from  the  QUIT  menu. 

This  takes  you  back  to  the  Xpmf  display.  Quit  from  it  simi¬ 
larly. 

6.  De-iconify  the  Satellite  Image. 

Remove  any  overlays  on  the  image  using  the  REFRESH 

menu. 

7.  Choose  Event  magnitude  data->Moment  Tensor  from  the 
OVERLAYS  menu. 

8.  Choose  Historical  magnitude  data->Moment  Tensor  from  the 
OVERLAYS  menu. 

9.  Choose  Surface  geology  from  the  OVERLAYS  menu. 

The  display  should  now  look  like  Figure  26. 

10.  Choose  from  the  QUIT  menu. 


1.  Choose  Yield  Estimation  from  the  main  FUNCTION  menu. 

2.  Click  on  MODEL  I  and  then  on  MODEL  3. 

The  display  should  look  like  Figure  27. 

3.  Set  the  weights  to  0.34, 0.33, 0.00  and  0.33,  respectively. 


4.  Click  on  MODEL  1  again. 

The  display  should  look  like  Figure  28. 

5.  Click  on  the  Save  button  in  the  lower  left  display. 

6.  Click  on  the  QUIT  button. 

1.  Choose  Statistical  Summary  from  the  main  FUNCTION  menu. 

2.  Click  on  Single  Event  Test  I  (AO, BO). 

After  a  few  seconds  an  icon  comes  up. 

3  Click  on  the  icon  to  see  the  test  results. 

4.  Click  on  the  QUIT  button. 
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Spreadsheet 


Quitting 
the  pro¬ 
gram 


Return  to  the  statistical  assessment  options  menu.  Use  steps 
2-4  to  run  Single  Event  Test  2  —  JVE  and  Single  Event  Test  3  --  Exchange 
Events. 

5.  Click  on  Single  Event  Test  1,2,3  Summary. 

6.  Click  on  the  icon  to  see  the  test  results. 

The  display  should  look  like  Figure  30. 

7.  Click  on  the  QUIT  button. 

Return  to  the  statistical  assessment  options  menu. 

8.  Set  AO  =  3.94  and  BO  =  0.81. 


Use  these  sliders  to  set  AO  and  BO 

9.  Click  on  Single  Event  Test  1  (A0,B0). 

When  the  icon  comes  up,  click  on  it  to  look  at  the  test 
results.  The  display  should  look  like  Figure  3 1 . 

10.  Click  on  the  QUIT  button  and  return  to  the  statistical  assess* 
ment  menu. 

11.  Click  on  the  QUIT  button  and  return  to  the  main  menu. 

1.  Choose  Spreadsheet  from  the  FUNCTION  menu. 

The  spreadsheet  shown  in  Figure  32  comes  up.  Use  the 
arrow  keys  to  move  around  the  spreadsheet. 

2.  Change  the  intercepts  to  4.47  (mi,),  4.47  (M]  g),  5.00  (Ms)  and 
4.48  (Mspec). 

Use  the  arrows  to  move  to  the  cell  and  type  (an  equal 
sign).  Then  type  the  number  and  hit  return. 

3.  Type  Q  to  quit  and  N  to  not  save  changes. 

Return  to  the  main  menu. 

1.  Choose  Quit  from  the  QUIT  menu. 
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Introduction 


This  section  provides  a  description  of  the  available  options  in  the 
YES  system.  It  describes  the  function  and  usage  of  every  module  and  its 
menus  and  buttons.  Following  is  a  review  of  basic  mouse  and  window 
techniques  required  by  the  user  to  work  with  YES. 

Mouse  Techniques 

Point  To  point,  position  the  cursor  on  something  you  want  to  choose. 

Click  To  click,  point  to  something  and  then  briefly  press  and  release  the 

left  mouse  button. 


Double  Click 


To  double-click,  point  to  something  and  then  rapidly  press  and 
release  the  left  mouse  button  twice. 


Drag  To  drag,  point  to  something,  press  and  hold  down  the  left  mouse 

button  as  you  move  the  cursor  to  another  position,  and  then  release  the 
mouse  button. 


Highlight 
a  button 


Choose  from 
a  pull-down 
menu 


To  highlight  or  turn  on  a  button,  position  the  pointer  on  the  button 
name  and  click  on  it. 

To  choo.se  an  option  from  a  pull-down  menu,  press  and  hold  down 
on  the  menu  button  name  to  pull  down  a  menu,  drag  the  highlight  to  the 
option  you  want  and  then  release  the  mouse  button. 
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Choose  from 
a  submenu 


Move  a 
window 


Iconify  a 
window 

De-lconify 
a  window 


To  choose  an  option  from  a  submenu,  first  pull  down  the  menu  that 
contains  the  submenu.  While  holding  down  the  mouse  button,  drag  diago¬ 
nally  to  the  submenu  indicator  (arrow)  to  pull  down  the  submenu.  Drag  the 
highlight  to  the  required  option  in  the  submenu  before  releasing  the  mouse 
button. 


Function  Selection 


Satellite  image 
World  Map 
Analyst  Station 
Magnitude  Mcasutcinent- 
Vield  Kstirnatioii: 
SiatiSiicd  ^oRiiuuy 
Spreadsfaibet 
Btertf  Summeary 


Submenu  indicator 


Window  Techniques 

To  move  a  window,  point  to  the  title  bar  of  the  window  to  be 
moved.  Click  and  hold  down  the  middle  mouse  button.  Drag  the  window  to 
the  desired  location  and  release  the  mouse  button. 

To  iconify  a  window,  point  to  the  small  box  in  the  far  left  of  the  top 
border  and  click  the  left  button.  The  iconified  window  can  be  moved  by 
clicking  and  holding  the  middle  button  inside  the  icon  and  dragging. 

To  de-iconify  a  window,  point  to  the  icon  and  click  on  the  left 
mouse  button. 


Tttle  bar  of  a  window 
Box  used  to  iconify  the  window 
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Main  Menu 


A 

B 


The  main  menu  provides  top  level  access  to  the  system.  It  consists 
of  six  menu  buttons  which  can  be  used  to  initiate  (SITE,  EVENT,  PHASE, 
FUNCTION)  or  terminate  (QUIT)  action  within  the  system,  or  to  view  on¬ 
line  information  regarding  the  operating  characteristics  and  parameters  of 
the  system  (HELP).  Choosing  any  of  these  buttons  causes  a  series  of  pull¬ 
down  menus  to  be  activated  from  which  the  various  system  options  can  be 
selected. 


A  -  the  main  menu  buttons 

B  -  the  information  bar  that  indicates  the  selected  site, 
event  and  phase 


Site 


The  SITE  menu  button  is  used  to  select  the  nuclear  test  site  at 
which  the  explosion  to  be  analyzed  was  detonated.  Currently,  the  choices 
are  limited  to  the  Soviet  Shagan  River  (Balapan)  and  Novaya  Zemlya  test 
sites.  Generally,  a  selection  must  first  be  made  from  this  menu  before  mak¬ 
ing  any  other  selections. 


Event 


Evta  Stie^i^Shacm 


The  EVENT  menu  button  is  used  to  .select  the  specific  explosion  to 
be  prtKCSsed  from  among  those  for  which  digital  seismic  data  are  currently 
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available  on  the  system. 

Phase 

The  PHASE  menu  button  is  used  to  select  from  among  the  six  seis¬ 
mic  arrivals  for  which  data  are  currently  available  on  the  system.  These 
consist  of: 

•  the  initial  P  and  Lg  arrivals  recorded  at  regional  distance  stations 
(Regional  P.  Lg), 

•  the  initial  short-period  teleseismic  P  arrivals  (Teleseismic  P), 

•  the  three  orthogonal  components  of  the  long-period  surface  waves 
(Rayleigh  ( vertical),  Rayleigh  ( radial),  and  Love). 

A  test  site  and  an  explosion  must  be  specified  from  the  SITE  and  EVENT 
menus  before  a  phase  is  selected. 

Function 

The  FUNCTION  menu  is  used  to  provide  access  to  the  eight  princi¬ 
pal  computational  and  analysis  modules  of  YES.  The  brief  overviews  of 
these  modules  which  are  presented  below  are  followed  in  subsequent  sec¬ 
tions  by  complete  descriptions  of  their  functionality.  Both  the  test  site  and 
explosion  of  interest  should  be  specified  before  selecting  from  this  menu. 

Satellite  Image 

This  module  provides  acce.ss  to  an  interactive  test  site  information 
interface  built  upon  displays  of  SPOT  satellite  images  of  the  Shagan  River 
and  Novaya  Zemlya  test  sites. 

World  Map 

This  module  provides  a  capability  to  display  the  map  locations  of 
stations  which  have  recorded  digital  seismic  data  from  the  selected  event. 


|Function|; 

RairtkWiScliectkm 
SktilUtte  Inn^ 
Woridklatp 
Amlyat  Station 
Aiagaitude  Measurement-^ 
ytel<l4*tlniaucih 
StatisiCjEeal  Summary 
Sprewbheet 
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Analyst  Station 

This  module  provides  the  capability  to  access  and  display  the  digital 
seismic  waveform  data  available  from  the  explosion  under  investigation  and  to 
interact  with  them  to  extract  the  parametric  data  required  for  determination  of 
the  various  magnitude  measures.  In  addition  to  the  test  site  and  event,  the  seis¬ 
mic  arrival  of  intere.st  must  be  selected  from  the  PHASE  menu  before  activat¬ 
ing  this  module. 

Magnitude  Measurement 

This  module  provides  acce.ss  to  the  algorithms  which  can  be  used  to 
semi-automatically  compute  network-averaged  magnitudes  from  the  wave¬ 
forms  and  parameters  corresponding  to  the  individual  seismic  arrivals  previ¬ 
ously  reviewed  by  the  operator  in  the  Analyst  Station  module.  At  the  present 
time,  five  .separate  magnitude  measures  ba.scd  on  tclesei.smic  P  waves  (m;,. 
Spectral  P),  long-period  Rayleigh  and  Love  waves  {Moment  Tensor),  and 
regional  P  and  Lg  waves  (Regional  P,  RMS  Lg)  can  be  estimated  using  this 
module. 

Yield  Estimation 

This  module  provides  access  to  a  statistical  module  which  permits  the 
operator  to  obtain  unified  yield  e.stimates  corresponding  to  the  multiple,  net¬ 
work-averaged  magnitudes  determined  using  the  Magnitude  Measurement 
module,  as  well  as  quantitative  measures  of  the  uncertainty  in  those  estimates. 

Statistical  Summary 

This  module  provides  access  to  a  variety  of  statistical  tests  which  can 
be  u.sed  to  quantitatively  assess  seismic  compliance  with  the  150  kt  threshold 
level  of  the  TTBT 
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Spreadsheet 

This  module  provides  the  analyst  with  the  capability  to  evaluate 
the  seismic  yield  estimate  for  the  current  event  in  a  standard  spreadsheet 
analysis  environment,  in  the  context  of  the  results  obtained  from  analyses 
of  previous  explosions  at  that  test  site,  as  well  as  any  available  independent 
yield  calibration  data  (c.g.,  CORRTEX). 

Event  Summary 

This  menu  item  provides  access  to  a  one  page  summary  of  the 
results  obtained  by  applying  the  YES  to  the  digital  seismic  data  recorded 
from  the  selected  explosion. 

Help 


In  this  and  all  subsequent  modules  to  be  described,  the  HELP  menu 
is  used  to  access  on-line  information  regarding  the  mode  of  activation  and 
functionality  of  the  various  menu  options  in  the  current  display. 


Quit 


In  this  and  all  subsequent  modules  to  be  described,  the  QUIT  menu 
is  used  to  terminate  the  operation  of  the  process  which  is  currently  runn  ng. 
Choosing  this  button  activates  a  pop-up  which  must  also  be  chosen  to  ver¬ 
ify  the  request  to  terminate  the  program. 
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Satellite  Image 

As  was  indicated  in  the  summary  overview  presented  in  the  preced¬ 
ing  section,  selection  of  the  Satellite  Image  option  from  the  main  FTJNC- 
TION  menu  provides  access  to  an  interactive  test  site  information  interface 
built  upon  displays  of  SPOT  satellite  images  of  the  Soviet  Shagan  River 
and  Novaya  Zemlya  test  sites.  For  the  purposes  of  this  application,  the 
image  for  the  Shagan  River  test  site  has  been  compressed  to  an  effective 
resolution  of  about  30  m  so  that  the  entire  test  site  can  be  viewed  on  a  sin¬ 
gle  screen.  Full  10  m  resolution  images  of  operator  designated  subareas  can 
be  accessed  using  the  Full  res  menu  option  described  below.  This  feature 
does  not  apply  to  Novaya  Zemlya,  since  only  20  m  resolution  SPOT  data 
are  currently  available  for  that  test  site.  In  either  case,  selection  of  Satellite 
Image  from  the  main  FUNCTION  menu  results  in  a  display  of  the  SPOT 
image  of  the  designated  test  site,  with  the  best  available  locations  of  previ¬ 
ous  explosions  at  that  site  shown  as  color-coded  square  overlays  and  with 
the  current  event  highlighted  by  a  yellow  diamond. 


A  -  satellite  image  menu  buttons 
B  -  information  bar 
C  -  side  menu 
O  -  current  event 


Reference 


Satellite  Image 


The  different  colors  are  used  to  differentiate  those  explosions  ?hout 
which  the  Soviets  have  published  data  in  the  open  literature  (blue)  from 
those  for  which  only  seismic  information  is  available  (red).  As  can  be  seen 
from  the  display,  selections  within  this  module  can  be  made  from  menu 
buttons  distributed  both  across  the  top  and  down  the  right  side  of  the  dis¬ 
play.  to  initiate  various  interactions  with  the  image  and  associated  data. 
These  will  now  be  described  in  order,  proceeding  from  left  to  right  and 
then  top  to  bottom. 


Images 


IMdtliaM  flf  aiiuvtitm 


This  menu  button  is  used  to  access  the  compressed  30  m  resolution 
SPOT  image  of  the  selected  test  site  {Satellite  image)  or  color-coded  repre¬ 
sentations  of  other  data  which  have  been  registered  to  the  SPOT  image.  For 
the  Shagan  River  test  site  these  include  the  DMA  topographic  data  (Sur¬ 


face  topography),  thickness  of  the  surface  alluvium  (Thickness  of  alluvium) 
and  the  depths  beneath  the  surface  to  the  various  subsurface  geologic  units 
which  have  been  mapped  (Depth  to  mesozoic.  Depth  to  paleozoic  and 
Depth  to  granite).  For  the  Novaya  Zemlya  test  site,  where  the  subsurface 
geology  is  less  well  known,  only  the  Surface  topography  (DMA)  image  is 
currently  available.  Where  appropriate,  color  scales  are  displayed  with  the 
images,  listing  numerical  values  of  the  displayed  variable  in  meters. 


Locations 


Lmoakn  tlwtfap 
Em*  iacallaM'(SK>1)' 
Evnat  localtim 

j22SESSL£2SS£^— 


The  LOCATIONS  menu  is  used  to  overlay  different  explosion  loca¬ 
tion  estimates  on  the  test  site  image  currently  being  displayed.  The  choice 
of  Event  locations(SPOT)  or  Event  locations(lSC)  from  this  menu  results 
in  the  overlay  of  the  corresponding  location  estimates  for  all  previous 
explosions  at  that  test  site.  The  choice  of  Compare  events  results  in  the 
overlay  of  the  best  available  locations  for  just  those  events  for  which  data 
are  available  for  a  particular  seismic  phase  and  station  selected  using  the 
Compare  function  from  the  Analyst  Station  module. 

The  location  symbols  can  be  used  to  obtain  information  about  the 
event  or  to  designate  the  event  for  further  proce.ssing.  For  example,  click- 
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ing  on  any  location  symbol,  causes  available  summary  information  about 
the  event  to  be  displayed  on  the  information  bar. 


Overlays 


QnwrUy  ifctecMwi 


CBnwt«««al 


evvnrausnirndtdaia 
metcrieat  m^itfudr  ihr  a  • 


DWBwd  oviBl»y» 


The  OVERLAYS  menu  is  used  to  access  and  view  a  variety  of  data 
which  have  been  registered  to  the  SPOT  satellite  images  of  the  two  test 
sites.  These  data  can  be  graphically  superimposed  on  any  of  the  available 
images  and  automatically  adjust  to  the  scale  of  the  selected  image. 


Current  event 


This  option  overlays  the  current  event  location  on  the  displayed 

image. 

Surface  geology 

This  option  overlays  a  graphical  representation  of  the  mapped  sur¬ 
face  geologic  features  for  the  selected  test  site  on  the  displayed  image  (see 
Figure  6  of  Chapter  2). 

Topography  contours 

The  Topography  contours  overlay  selection  can  be  used  to  superim¬ 
pose  a  contour  representation  of  the  DMA  topographic  data  for  the  selected 
test  site  on  either  the  compressed  or  full  resolution  SPOT  images  or  on  the 
color-coded  surface  topography  representation  from  the  IMAGES  menu. 
The  contours  are  labeled  with  the  corresponding  values  of  elevation  above 
sea  level  in  meters  (see  Figure  7  of  Chapter  2). 


Event  magnitude  data 


Event  magnltutie  data 
Mbiesitinal 
.Mb-MLg 
Mameit  tensor 
Mb-Mo _ _ 


Selection  of  the  Event  magnitude  data  option  generates  a  pull 
down  menu  listing  the  various  magnitude  related  overlays  available  for  the 
current  event.  Magnitudes  must  have  already  been  estimated  for  this  event 
before  these  overlays  are  acce.ssed.  If  the  required  magnitudes  are  not 
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available  for  a  given  selection,  a  message  will  appear  on  the  screen  stating 
“overlay  file  not  available.”  The  Mb  residual  option  produces  an  overlay  to 
the  displayed  image  listing  the  difference  between  the  current  event  net¬ 
work-averaged  mj,  value  and  the  mj,  value  obtained  for  that  event  at  an 
individual  station  specified  in  the  Magnitude  Measurement  module,  plotted 
at  the  current  event  location.  Similarly,  selection  of  the  Mb  -  MLg  or  Mb  - 
Mo  option  produces  overlays  listing  the  difference  between  the  current 
event  network-averaged  m^,  value  and  a  corresponding  single  station  Ml^ 
value  (  see  Figure  20  of  Chapter  2)  designated  in  the  Magnitude  Measure¬ 
ment  module,  or  network-averaged  surface  wave  isoU'Opic  moment  esti¬ 
mate  (Mq),  respectively.  Selection  of  the  Moment  tensor  option  from  this 
menu  produces  an  overlay  of  a  graphical  representation  of  the  surface  wave 
moment  tenor  inversion  solution  at  the  current  event  location  (see  Figure 
26  of  Chapter  2).  This  display  consists  of  two  concentric  circles  enclosing  a 
straight  line  segment,  where  the  ratio  of  the  diameters  of  the  circles  is  pro¬ 
portional  to  the  ratio  of  the  isotropic  (light  blue)  and  tectonic  (yellow) 
moments  and  the  yellow  line  segment  is  oriented  parallel  to  the  inferred 
strike  of  the  tectonic  release  component. 


Historical  magnitude  data 


Mbmidaato  -• 

Mb-MI«(NORSAIt) 

MttncKtvaMT 


Selection  of  Historical  magnitude  data  from  the  OVERLAYS 
menu  activates  a  pull  down  menu  which  provides  access  to  overlay  repre¬ 
sentations  of  previous  experience  with  the  various  magnitude  related 
parameters  corresponding  to  the  current  Event  magnitude  data  selections 
which  were  just  described.  In  this  case,  selection  of  Mb  residuals  activates 
a  secondary  pull  down  menu  listing  the  various  USAEDS  stations  for 
which  the  previously  observed  differences  between  network-averaged  and 
single  station  m^  values  have  been  contoured  as  a  function  of  event  loca¬ 
tion  within  the  test  site  (currently  only  available  for  the  Shagan  River  test 
site).  Similarly,  selection  of  Mb  -  MLg  (NORSAR)  produces  an  image  over¬ 
lay  of  the  differences  between  the  network-averaged  mj,  values  and  corre¬ 
sponding  NORSAR  Lg  values  in  either  contour  (Shagan  River)  or  discrete 
event  (Novaya  2^mlya)  form.  The  Moment  tensor  option  produces  an  over¬ 
lay  of  graphical  rcpre.scntations  of  all  previous  .surface  wave  moment  ten¬ 
sor  solutions  for  explosions  at  the  selected  test  site,  where  the  i.sotropic  and 
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tectonic  components  are  represented  by  dark  blue  and  red  circles,  respec¬ 
tively,  tf)  differentiate  them  from  the  corresponding  current  event  solution 
(see  Figure  26  of  Chapter  2).  Selection  of  Mh  -  Mo  produces  an  overlay  of 
the  differences  between  the  network-averaged  mj,  and  values  for  all 
previ(ms  events  for  which  both  these  magnitude  measures  are  available. 


Subsurface  contours 


fKimmtact  eaitwn 
TMrlmwtciWtoaw 
Dcptb  to  Up  «(  Memglc 
Di!|Mi  to  of  paleoKctc 

nopH>wmip«fira«tu 


The  Subsurface  contours  option  in  the  OVERLAYS  menu  provides 
the  capability  to  overlay  contours  of  the  various  subsurface  geologic  layer 
thicknesses  and  depths  on  a  displayed  image,  generally  the  corresponding 
image  from  the  IMAGES  menu.  The  contours  are  labelled  with  the  corre¬ 
sponding  numerical  values  of  the  displayed  variable  in  meters. 


Location  differences 


LocoUtPdmtettBOM 

PnfemdtoSPOT 

PicfencdtoEC 

PrtfcntdtoBocharov 


The  Location  differences  option  provides  access  to  a  set  of  misloca- 
tion  vector  overlays  consisting  of  directed  line  segments  pointing  from  the 
preferred  to  an  operator-specified  set  (e.g.,  SPOT,  ISQ  of  explosion  loca¬ 


tion  estimates.  The  preferred  locations  are  those  which  are  thought  to  be 
most  accurate  from  among  those  available,  which  in  the  classified  environ¬ 


ment  are  generally  the  NPIC  photo  locations. 


Digitized  overiays 


DIgWzpd  owatoyt 
Corrent  overlay 
Shagan  wMlvIafaoa 


The  Digitized  overlays  option  provides  access  to  an  extensible  list 
of  overlay  files  which  have  been  generated  by  the  operator  using  the  Digi¬ 
tize  feature  to  be  described  below.  If  a  feature  has  been  digitized  from  the 
SPOT  image  in  the  current  session,  it  can  be  overlaid  on  the  displayed 
image  by  selecting  Current  overlay  from  this  menu.  The  results  of  any  dig¬ 
itizations  which  have  been  saved  to  a  permanent  file  by  the  operator  in  pre¬ 
vious  sessions  (e.g.,  Shagan  subdivisions)  can  also  be  accessed  from  this 
menu  and  overlaid  on  any  of  the  available  images. 
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Refresh 


Benov«0vecl^  This  menu  option  is  used  to  remove  any  overlays  {Remove  over- 

itemove  tocattotts  I 

lays)  or  location  symbols  {Remove  locations)  which  are  currently  superim¬ 
posed  on  the  displayed  image. 

Function 

luvGJae  a  fUnctloaT 

Anaiy«*  rfaHmt  This  menu  button  is  used  to  initiate  processes  and  create  permanent 

S^  ovgriay  I  overlay  files. 

Analyst  station 

Selection  of  Analyst  station  initiates  a  process  by  which  seismic 
waveforms  recorded  at  a  specified  station  from  events  selected  using  the 
side  menu  Compare  function  described  below  are  retrieved  from  the  data¬ 
base  and  displayed  on  the  screen  in  the  standard  Analyst  station  format  (cf. 
Analyst  Station  description). 

Save  overlay 


Selection  of  Save  overlay  initiates  a  process  by  which  a  digitized 
overlay  can  be  saved  to  a  permanent  file  to  appear  on  the  Digitized  overlays 
menu  on  all  sub.sequeni  restarts  of  the  system.  This  feature  is  used  in  con¬ 
junction  with  the  Digitize  button  in  the  side  menu.  It  is  implemented 
through  a  pop-up  menu  shown  below  which  allows  the  operator  to  type  in  a 
name  for  the  file  to  be  saved  (e.g.,  Shagan  subdivisions)  and  then  save  it  to 
a  permanent  file  by  clicking  on  the  Save  button.  If  at  any  time  during  the 
process  the  operator  decides  against  saving  the  file,  the  sequence  can  be 
terminated  by  clicking  on  the  Cancel  button. 


Save 


Cancel 


Type  In  overlay  name  (<  25  charactersM 


A  -  text  window  used  to  type  in  a  name. 
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I  Ut/lon  t 
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t  guUirw  t 
lCro»»«>gt 


TCcrofO 


Side  Menu 

The  side  menu  items  on  the  satellite  image  display  provide  addi¬ 
tional  options  which  permit  the  operator  to  interact  with  the  display  in  a 
variety  of  different  ways.  A  notable  feature  of  these  buttons  is  that  they  ini¬ 
tiate  processes  which  remain  active  while  that  button  is  turned  on.  Conse¬ 
quently,  they  must  be  turned  off  when  the  associated  processing  has  been 
completed. 

Lat/lon 

When  this  button  is  highlighted,  clicking  on  any  point  on  the  image 
causes  the  latitude  and  longitude  of  that  point  to  be  displayed  in  the  infor¬ 
mation  bar. 

Digitize 

When  the  Digitize  button  is  highlighted,  clicking  on  any  point  or  a 
series  of  points  on  the  image  causes  the  latitude  and  longitude  of  the 
point(s)  to  be  added  to  a  temporary  file.  When  digitization  of  the  selected 
feature  has  been  completed,  the  file  is  closed  by  clicking  on  the  Digitize 
button  again.  Any  subsequent  selection  of  Current  overlay  from  the  Digi¬ 
tized  overlays  in  the  OVERLAYS  menu  will  cause  straight  line  segments  to 
be  drawn  between  these  points  in  the  order  in  which  they  were  selected  and 
overlaid  on  the  image.  If  desired,  this  information  can  also  be  .saved  to  a 
permanent  named  file  using  the  Save  overlay  option  from  the  Satellite 
Image  FUNCTION  menu. 

Fuii  res 

When  the  Full  res  button  is  highlighted,  in  conjunction  with  the 
Shagan  River  test  site  display,  clicking  on  a  point  on  the  image  causes  the 
full  10  m  resolution  SPOT  image  data  to  be  di.splayed  for  the  area  sur¬ 
rounding  that  point.  Once  a  full  re.solution  image  has  been  di.splayed,  full 
resolution  images  of  adjacent  regions  can  be  acces.sed  by  pointing  the  cur¬ 
sor  near  the  closest  border  of  the  display  and  clicking  on  the  mouse  button. 
To  access  the  full  re.solution  image  for  a  new  area  not  adjacent  to  the  cur- 
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rent  display,  select  Satellite  image  from  the  IMAGES  menu  and  select  the 
new  location.  Upon  completion  of  the  review  of  the  full  resolution  data, 
turn  off  the  Full  res  button  and  proceed. 

Cross  sec 

When  this  button  is  on,  a  vertical  subsurface  geologic  section  along 
a  line  between  any  two  points  on  the  image  can  be  created  and  displayed  by 
sequentially  selecting  the  points  with  the  mouse.  The  Quit  button  on  the 
cross  section  display  is  used  to  refresh  the  underlying  image  display,  from 
which  another  cross  section  can  be  initiated,  or,  upon  completion,  the  pro¬ 
cess  can  be  terminated  by  turning  off  the  Cross  sec  button. 

Compare 

When  this  button  is  on,  events  can  be  selected  with  the  mouse  from 
among  those  displayed  on  the  image  using  the  Compare  events  option  from 
the  LOCATIONS  menu.  Clicking  on  any  location  symbol  causes  the  event 
location  to  be  highlighted  by  a  diamond  and  available  summary  informa¬ 
tion  about  the  event  to  be  extracted  from  the  database  and  displayed  on  the 
information  bar.  If,  on  the  basis  of  this  information,  the  operator  decides  to 
select  this  event,  a  second  click  of  the  mouse  button  will  cause  that  event 
identification  label  to  be  added  to  a  temporary  file  and  the  highlight  dia¬ 
mond  will  turn  white.  Additional  events  can  be  selected  in  the  same  fashion 
and,  upon  completion,  the  file  is  closed  by  turning  off  the  Compare  button. 
Subsequent  selection  of  Analyst  station  from  the  FUNCTIONS  menu 
causes  the  designated  station  waveforms  corresponding  to  the  selected 
events  in  this  file  to  be  retrieved  from  the  databa.se  and  displayed  on  the 
screen,  with  the  current  event  recording  at  the  top  and  the  selected  event 
waveforms  below  it  in  the  order  in  which  they  were  selected. 

Contrast  and  Brightness 

The  slider  bars  at  the  bottom  of  the  side  menu  can  be  used  to  inter¬ 
actively  adjust  the  contra.st  (C)  and  brightness  (B)  of  the  displayed  SPOT 
image.  Clicking  on  the  C  or  B  buttons  returns  the  display  to  the  default  set¬ 
ting  for  that  parameter. 
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World  Map 

This  module  provides  a  capability  to  display  the  map  locations  of 
stations  which  have  recorded  digital  seismic  data  for  the  selected  event. 
The  maps  are  plotted  as  azimuthal  equidistant  projections  of  the  globe  (A  < 
100°)  centered  on  the  selected  test  site. 


A 

B 


A  -  world  map  menu  buttons 


B  -  information  bar 


C  -  current  test  site 


Phase 


This  menu  button  is  used  to  select  the  seismic  phase  for  which  the 
database  will  be  searched  for  stations  which  recorded  data  from  the  cur¬ 
rent  ’vent.  The  locations  of  these  stations  are  then  automatically  plotted 
on  the  map  display  together  with  the  standard  three  or  four  letter  station 
abbreviations. 


61 


Reference 


World  Map 


Station 


1  sumdim  1 

[  ManailF  i  l«  1 

lUylCl(fa<ndU> 

HlA 

IaIhd 

IGAX 

MV 

lOBK 

KOUO 

WMQ 

This  menu  button  is  used  to  select  a  station  from  among  those 
which  have  recorded  data  for  the  designated  phase  and  event  Selection  of  a 
station  causes  a  straight  line  to  be  drawn  on  the  map  between  the  test  site 
and  the  station  which,  in  this  projection,  corresponds  to  the  great  circle 
path  between  these  two  points. 


Refresh 


This  menu  button  is  used  to  remove  the  station  location  overlays 
from  the  map  display  prior  to  selection  of  another  phase. 


62 


Reference 


Analyst  Station 


Analyst  Station 

Selection  of  Analyst  Station  from  the  main  FUNCTION  menu  fol¬ 
lowing  the  selection  of  an  event  and  seismic  phase  produces  a  display  of 
the  corresponding  digital  seismic  waveform  data  in  a  format  which  permits 
the  operator  to  interact  with  them  to  extract  the  parametric  data  required  for 
determination  of  the  various  magnitude  measures  for  the  current  event.  The 
default  setting  provides  for  the  waveforms  to  be  di.splayed  in  order  of 
increasing  station  epicentral  distance,  with  six  waveforms  per  screen  or 
“page.”  The  mode  of  di.splay  varies  with  the  selected  phase,  with  teleseis- 
mic  and  regional  P  segments  denoted  (P  and  Pn)  displayed  as  a  function  of 
time  (from  10  seconds  before  to  15  seconds  after  signal  onset),  while  the  Lg 
and  surface  wave  segments  are  displayed  as  a  function  of  group  velocity 
(10  to  2.5  km/sec  for  Lg;  5  to  2.5  km/sec  for  the  Rayleigh  and  Love 
phases). 


A  -  analyst  station  menu  buttons  C  -  waveform  information 

B  -  the  information  bar  D  -  arrival  pick 


As  is  illustrated  in  the  display,  each  waveform  segment  is  annotated 
to  indicate  the  event,  station,  phase,  epicentral  distance,  approximate  single 
station  mj,  value  and  a  waveform  quality  assignment  which  can  be  interac¬ 
tively  set  using  the  TOOLS  menu  function  Set  quality  to  be  described 
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below.  The  vertical  line  segments  on  these  displays  marks  the  estimated 
signal  onset  time  which  also  can  be  changed  interactively  using  the 
TOOLS  menu  function  Pick  arrival  to  be  described  below. 

Next 

This  menu  button  is  used  to  page  forward  through  the  waveform 
displays,  six  waveforms  at  a  time,  in  order  of  increasing  station  epicentral 
distance. 

Previous 

This  menu  button  is  used  to  page  backwards  from  the  current  dis¬ 
play  through  the  available  waveforms  from  stations  at  smaller  epicentral 
distances. 

Tools 

This  menu  button  is  used  to  access  the  various  signal  processing 
and  parameter  extraction  functions  which  have  been  provided  to  facilitate 
the  evaluation  of  the  displayed  seismic  data  and  the  determination  of  the 
various  seismic  magnitude  measures.  Once  a  selection  has  been  made  from 
this  menu,  the  corresponding  function  remains  active  until  another  choice 
is  made  or  until  a  new  analyst  station  is  initiated. 

Pick  arrival 

This  tool  permits  the  operator  to  change  the  estimate  of  signal  onset 
time  by  pointing  at  the  location  of  the  new  estimate  on  that  trace  and  click¬ 
ing  once.  Prior  to  analyst  review,  the  default  onset  times  are  the  predicted 
Herrin  68  travel  times  for  P  and  Pp  and  a  group  velocity  of  about  3.5  km/ 
sec  for  the  Lg  and  surface  wave  phases.  Upon  completion  of  this  review, 
the  revised  signal  onset  estimates  can  be  saved  to  the  database  using  the 
Save  feature  from  the  FUNCTION  menu  described  below.  Subsequent  ana¬ 
lyst  station  displays  of  these  data  will  then  show  the  revised  arrival  esti¬ 
mates. 


To>o| 

PIde  arrival 
Setqtuitcy 
Deniettt/dctiaal 
Fllt«r 

Neiwr  tnatniBiiCaa  - 

Expaul 

mb 

Ccnqimre 
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DataquaUty 
Excellettf 
Good 
VJsiHe 
bivlsibie 
Broken 
bdeifeiing 
Qipped _ 


Set  quality 

Selection  of  this  tool  activates  a  pull  down  menu  listing  a  variety  of 
quality  categories  which  can  be  assigned  to  the  displayed  data.  If  no  prior 
analyst  assignment  has  been  saved  to  the  database,  the  trace  will  be  dis- 
plascd  initially  with  a  Not  Checked  label.  Assignment  of  a  revised  quality 
can  be  accomplished  by  selecting  a  category  from  this  menu  and  then 
clicking  on  any  trace  to  which  that  category  applies.  Category  assessments 
are  subject  to  operator  judgment,  but  generally  Excellent  is  associated  with 
apparent  signal-to-noisc  ratios  of  greater  than  about  five.  Good  with  signal- 
lo-noise  ratios  in  the  range  from  about  two  to  five  and  Visible  to  identified 
.signals  with  signal-to-noi.se  ratios  of  less  than  about  two.  The  la.st  four  cat¬ 
egories  {Invisible.  Broken.  Interfering.  Clipped)  are  generally  used  to  des¬ 
ignate  data  which  will  not  be  processed  to  obtain  magnitude  estimates. 
Invisible  signities  that  the  displayed  dam  contains  no  evident  signal,  while 
Broken  indicates  an  in.strument  malfunction  and  Interfering  denotes  con¬ 
tamination  of  the  displayed  data  by  arrivals  from  a  different  event.  The  cat¬ 
egory  Clipped  is  assigned  only  if  that  portion  of  the  displayed  data 
containing  the  phase  arrivals  to  be  used  in  the  corresponding  magnitude 
determination  shows  evidence  of  clipping.  New  quality  assignments  can  be 
saved  to  the  databa.se  at  any  time  by  using  the  Save  feature  from  the 
FUNCTION  menu  de.scribcd  below. 

Demean/detrend 

This  tool  permits  the  operator  to  remove  any  DC  bias  or  linear  trend 
from  traces  which  are  subsequently  selected  with  the  mouse.  The  resulting 
di.splayed  time  series  have  mean  amplitudes  of  zero. 

Filter 


This  tool  permits  the  operator  to  bandpass  filter  (Butterworth)  any 
subsequently  designated  trace  between  interactively  set  low  frequency  and 
high  frequency  limits.  These  frequency  bounds  are  set  with  slider  bars 
.icces.od  from  a  pop-up  window  (.see  Figure  14  of  Chapter  2),  as  are  the 
4opes  (i.e.,  pole  numbers)  of  the  respon.se  roll-off  outside  the  specified 
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band.  When  the  bandpass  filter  processing  has  been  completed,  the  interac¬ 
tive  slider  window  can  be  cancelled  by  selecting  the  Exit  button. 

New  Instrument 


New  fiiatniineM 
ww«ai8i> 
wwtMEaliP 
sro^ 


This  tool  permits  the  operator  to  analytically  transform  any  subse¬ 
quently  designated  trace  into  the  trace  which  would  have  been  observed  at 
that  station  if  the  data  had  been  recorded  on  any  one  of  three  standard 
instruments  (i.e.,  WWSSN,  SRO,  USAEDS)  whose  response  characteris¬ 
tics  are  documented  in  the  EFCRES  instrument  response  file  at  the  DARPA 
Center  for  Seismic  Studies. 


Expand 


This  tool  permits  the  operator  to  produce  a  full  page  display  of  any 
subsequently  designated  trace.  For  P  and  Pp,  both  the  amplitude  and  time 
scales  are  expanded,  with  the  displayed  time  segment  extending  from  four 
seconds  before  to  six  seconds  after  the  estimated  signal  onset  time.  For  the 
Lg  and  surface  wave  displays,  only  the  amplitude  .scale  is  expanded.  This 
feature  is  most  commonly  used  in  conjunction  with  the  tool  to  facilitate 
the  determination  of  peak-to-peak  amplitude  and  period  of  the  pulse  used 
to  define  the  single  station  P  wave  magnitude  value.  However,  any  of  the 
other  TOOLS  menu  functions  can  be  accessed  from  the  expanded  display 
and  applied  to  the  data.  When  processing  of  the  expanded  data  has  been 
completed,  the  display  can  be  terminated  by  selecting  the  Exit  button. 

mt, 


This  tool  allows  the  operator  to  use  the  mnu.se  to  draw  a  rectangle 
on  a  selected  trace  so  as  to  define  the  peak-to-peak  amplitude  and  half 
period  of  the  selected  cycle  of  initial  P  or  P^  motion  to  be  used  in  the  esti¬ 
mation  of  the  corresponding  single  station  magnitude  measure.  This  is 
accomplished  by  clicking  the  cursor  at  the  peak,  holding  down  the  mouse 
and  dragging  the  cursor  to  the  trough  of  the  selected  cycle.  Upon  release  of 
the  button  the  approximate  peak-to-peak  ground  motion  amplitude  and 
a.s.sociated  period  of  the  pulse  is  automatically  computed  and  displayed  on 
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the  inlbrmation  line  ol  the  analyst  station  display.  After  this  process  has 
been  repeated  for  each  usable  trace,  the  amplitudes  and  periods  can  be 
saved  to  a  file  for  subsequent  input  to  the  corresponding  magnitude  estima¬ 
tion  module  using  the  Interface  file  feature  from  the  FUNCTION  menu 
described  below. 

Compare 

This  tool  allows  the  operator  to  select  a  trace  and  initiate  a  process 
in  which  the  on-line  seismic  database  is  searched  for  all  other  events  at  the 
selected  test  site  for  which  data  for  the  displayed  phase  were  recorded  at 
the  designated  station.  A  file  identifying  the  events  meeting  these  criteria  is 
then  automatically  created  and  used  U'  update  the  file  u.sed  by  Compare 
events  in  the  LOCATIONS  menu  of  the  Satellite  Image  module.  The  map 
locations  of  the.se  events  can  then  be  overlaid  on  the  SPOT  satellite  image 
of  the  test  site  and  the  designated  station  data  corresponding  to  selected 
events  can  be  displayed  and  compared  to  the  data  recorded  at  that  station 
from  the  current  event  using  the  Compare  and  Analyst  station  functions 
from  the  Satellite  Image  menus. 

Clone 


This  tool  allows  the  operator  to  select  any  displayed  trace  and  cre¬ 
ate  a  color-coded  (red)  copy  of  that  trace  which  can  be  dragged  with  the 
mouse  to  any  position  on  the  screen  to  permit  detailed  comparison  with 
other  displayed  waveforms  (see  Figure  16  of  Chapter  2).  This  process  is 
applicable  to  any  P  or  Pp  display  or  to  any  analyst  station  display  of  com¬ 
pare  event  data  recorded  at  a  fixed  station.  However,  it  should  be  noted  that 
Lg  or  surface  wave  recordings  from  stations  at  different  distances  cannot  be 
directly  compared  in  this  fashion  because  they  are  di.splayed  as  a  function 
of  group  velocity  rather  than  time. 
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Function 

This  menu  button  is  used  to  save  processing  results  to  tempo. ;i:y 
(Interface  file)  or  permanent  (Save)  files,  or  to  undo  (Restore)  the  results  of 
any  signal  poKessing  which  has  been  applied  to  the  displayed  waveforms. 


Interface  file 


The  Interface  file  option  enables  the  operator  to  save  any  interac¬ 
tively  determined  parametric  data  to  temporary  files  which  will  subse¬ 
quently  be  used  as  input  to  the  appropriate  magnitude  determination 
modules.  More  specifically,  once  the  amplitudes  and  periods  for  a  particu¬ 
lar  P  phase  (P  or  P^)  have  been  determined  for  each  usable  recording  using 
the  nif,  function  from  the  TOOLS  menu,  the  results  can  be  saved  for  further 
processing  by  invoking  the  Interface  file  function. 

Save 


The  Save  option  enables  the  operator  to  save  any  revised  determi¬ 
nations  of  data  quality  (Set  quality)  and  signal  onset  times  (Pick  arrival)  to 
the  permanent  database.  Once  this  function  has  been  invoked,  any  subse¬ 
quent  displays  of  these  same  data  in  the  analyst  station  format  will  incorpo¬ 
rate  the  revised  estimates  of  these  parameters. 

Restore 

The  Restore  option  allows  the  operator  to  remove  the  effects  of  any 
liltering  or  instrument  translormations  conducted  during  the  current  ses¬ 
sion  to  return  the  displayed  waveforms  to  their  initial  form.  However,  any 
revisions  to  data  quality  and  signal  onset  estimates  made  during  the  current 
session  are  retained  during  this  process  and,  consequently,  the  Restore 
option  can  also  be  used  to  realign  the  displayed  traces  to  coincide  with  any 
revised  signal  onset  time  markers. 
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Reference  Ma  gnitude  Measurement 

Magnitude  Measurement 

This  module  provides  the  capability  to  estimate  network-averaged 
seismic  magnitude  measures  from  the  waveform  and  parameter  data  corre¬ 
sponding  to  the  individual  seismic  arrivals  previously  reviewed  by  the 
operator  using  the  analyst  station  module.  At  the  present  time,  five  separate 
magnitudes  based  on  regional  P  and  Lg  waves  {Regional  P,  RMS  Lg), 
lelescismic  P  waves  {mh.  Spectral  P)  and  long-period  Rayleigh  and  Love 
waves  {Moment  Tensor)  can  be  estimated  using  this  module. 

Regional  P 


Regional  P  Magnitude 

MlPa)  -  C.025 


A  -  network-averaged  magnitude 
B  -  standard  deviation 
C  -  mean  magnitude  estimate 

Used  to  estimate  single  station  and  network-averaged  regional  P 
phase  magnitudes  (  MiPp) )  from  the  amplitudes  and  periods  which  were 
saved  in  the  corresponding  interface  file  during  the  analyst  station  review 
of  the  regional  P  phase  data.  The  resulting  display  lists  the  network-aver¬ 
aged  magnitude  value  (  M(Pp) )  and  the  as.sociated  standard  deviation  of 
the  data  (o)  and  shows  the  corresponding  single  station  magnitude  esti¬ 
mates  (filled  circles)  with  respect  to  the  mean  and  plus  and  minus  la  lev- 
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Reference 


Magnitude  Measurement 


els  (horizontal  lines).  At  the  same  time,  the  (  M(Pn) )  value  is  automatically 
written  to  a  temporary  file  and  is  designated  as  the  regional  P  wav.  magni¬ 
tude  for  the  current  event.  The  menu  options  on  this  display  can  bo  u'  .o,  to 
interact  with  and  revise  this  fully  automatic  solution. 

Add/Delete 

Selection  of  the  Add/Delete  option  activates  a  pulldown  menu  list¬ 
ing  the  individual  stations  contributing  to  the  network-averaged  magnitude 
estimate.  If  the  operator  wishes  to  eliminate  a  station  from  the  average 
(e.g.,  because  of  large  deviation  from  the  mean  \  ue),  highlighting  that 
station  name  in  the  menu  causes  that  station  to  be  .  .linated  from  future 
estimates  of  the  network-averaged  value.  This  process  can  be  repeated  to 
delete  any  number  of  additional  stations. 

Recalc 

Once  a  desired  subset  of  stations  has  been  defined  through  this  pro¬ 
cess,  the  Recalc  menu  option  can  be  selected  to  automatically  generate  a 
revised  magnitude  estimate  and  corresponding  display.  At  any  point  in  this 
review,  an  eliminated  station  can  be  added  back  into  the  averaging  process 
by  re-selecting  it  from  the  Add/Delete  station  list,  at  which  point  the  high¬ 
light  on  that  station  name  will  be  eliminated. 

Save 


If  a  preferred  value  of  M(P„)  is  identified  during  this  review  pro¬ 
cess  it  can  be  saved  to  the  current  event  magnitude  file  using  the  Save  menu 
option,  at  which  point  it  replaces  the  complete  network  value  which  was 
automatically  written  to  this  file  after  the  initial  calculation. 

Residual 

.Selection  of  the  Residual  option  activates  another  pull  down  menu 
with  the  same  list  of  stations,  any  one  of  which  can  be  selected.  This  ini¬ 
tiates  a  prcKoss  by  which  the  difference  between  the  current  event  network- 


Me  ;nitude  Measurement 


averagcil  mj,  value  i,  rn(,)  and  the  uesignated  station  M(Pn)  value  is  com¬ 
puted  and  added  to  the  Event  magnitude  data  overlay  menu  in  the  Satellite 
Image  module,  where  it  can  be  superimposed  on  the  image  at  the  current 
event  location  and  compared  with  values  of  this  difference  obtained  for 
previous  events  at  that  test  site  (if  available).  Note  that  an  niii,  value  must 
have  already  been  estimated  for  the  current  event  before  this  option  can  be 
exercised. 

RMSLg 

This  module  is  used  to  estimate  single  station  and  network-aver¬ 
aged  regional  Lg  phase  magnitudes  (M(Lg))  from  the  vertical  component 
Lg  signals  recorded  at  stations  for  which  the  data  were  found  to  be  usable 
in  the  analyst  station  review  process.  In  computing  this  magnitude  mea¬ 
sure,  the  data  are  first  bandpass  filtered  in  a  band  which  typically  extends 
from  about  0.6  to  3.0  Hz  and  then  the  peak  of  the  RMS  amplitude  level  of 
the  filtered  signal  is  computed  over  a  group  velocity  interval  extending 
from  about  3.6  to  3.0  km/sec.  For  array  stations  such  as  NORSAR,  the 
M(Lg)  value  is  determined  as  a  logarithmic  average  of  the  peak  RMS 
amplitude  values  at  the  various  array  sites.  Since  this  is  a  computationally 
intensive  process,  it  is  generally  run  in  the  background  while  the  operator 
conducts  other  tasks.  When  the  process  has  been  completed,  an  RMS  Lg 
icon  will  appear  on  the  screen,  where  it  can  be  activated  by  the  operator  at 
his  discretion.  The  resulting  magnitude  display  and  associated  functionality 
RMS  Lg  icon  are  identical  to  those  described  above  for  the  Regional  P  phase  module 
with  M(Lg)  replacing  M(Pn).  As  with  the  Regional  P  phase  module,  an 
mj,  value  must  have  already  been  estimated  for  the  current  event  before 
the  RMS  Lg  Residual  menu  option  is  exercised. 

mtj 


Reference 


This  module  is  used  to  estimate  single  station  and  network-aver¬ 
aged  teleseismic  P  wave  magnitudes  from  the  amplitudes  and  periods 
which  were  saved  in  the  corresponding  interface  file  during  the  analyst  sta¬ 
tion  review  of  the  tele.seismic  P  phase  data.  The  resulting  magnitude  dis- 
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Magnitude  Measurement 


play  and  associated  functionality  are  identical  to  those  described  above  for 
the  Regional  P  phase  module  with  replacing  MCPp).  Thus,  in  opli- 
caticn  to  selection  of  a  station  from  the  j’rsidual  menu  initiates  a  nro- 
cess  by  which  the  difference  between  the  current  event  network-averaged 
mjj  value,  rn^,,  and  the  designated  station  m^  value  is  computed  and  added 
to  the  Event  magnitude  data  overlay  menu  m  the  Satellite  Image  module. 
At  the  present  time,  this  value  can  be  compared  to  previous  experience 
only  for  US  AEDS  station  recordings  of  Shagan  River  explosions  for  which 
these  single  station  m^,  residuals  have  been  contoured  as  a  function  of  event 
location  within  the  test  site. 

Spectral  P 

This  module  is  used  to  estimate  network-averaged  teleseismic  P 
wave  spectra  and  associated  spectral  P  magnitudes  (M(Pspec))  and  yield 
estimates  through  processing  of  teleseismic  P  phase  data  recorded  at  sta¬ 
tions  for  which  frequency  dependent  station  correction  factors  have  previ¬ 
ously  been  estimated  through  off-line  analysis  Once  Spectral  P  has  been 
selected  from  the  Magnitude  Measurement  menu,  a  separate  process  is  ini- 


A  -  Spectral  P  menu  buttons 
B  -  Stations  'w  indow 
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Mc  ^nitude  Measurement 


tiaied  with  its  own  display  and  associated  menu.  In  this  initial  display,  the 
usalih  stations  for  which  leleseismic  P  phase  data  are  available  are  listed  in 
the  Stations  window  in  the  lower  left  hand  corner  of  the  display,  together 
with  A(hl  and  Delete  buttons  which  can  be  used  by  the  operator  to  add  or 
delete  stations  from  the  proce.ssing  sU'cam  at  any  stage.  The  proce.ssing  is 
accomplished  by  accessing  the  menu  options  in  the  following  .sequence. 


NBF 


A  -  P  wave  signal 
B  -  signal  spectrum 
C  -  noise  spectrum 


This  module  permits  the  operator  to  estimate  and  display  single  sta¬ 
tion  P  wave  signal  and  noi.se  .spectra  through  narrowband  filtering  of  the 
waveform  data  over  the  band  from  0.50  to  2.25  Hz.  If  upon  review  of  these 
spectra  the  operator  concludes  that  the  data  from  a  particular  station  are 
unu.sable,  that  station  can  be  deleted  from  further  proce.ssing  by  .sequen¬ 
tially  .selecting  the  Delete  and  station  buttons  from  the  Stations  window. 
( )nce  the  review  of  the  .spectra  fr.r  a  given  station  is  completed,  .selection  of 
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Magnitude  Measurement 


Coniiriic  t'ioin  lUcaj  ;jjIi  i  :  in  ''!v>  i.'u;  scicen  and  brings  up 

the  special  dispk.v  r-nr.  .'.p-  i>n  ;  '  .  hr  rx.  xiaiion.  When  Co-  t'nue  is 
selC'j'ea  Icliownig  ii;c  .asi  Maaiai,  I'v  :■  tr  ;  ..  lion  is  automaiicai!>  lercni- 
nateu.  signaling  the  end  of  the  single  siaiioii  processing. 

Avgspec 


Seleciion  ol  liifs  menu  opiion  initiates  a  process  by  which  fre¬ 
quency  dependent  station  correction  factors  are  applied  to  the  single  station 
spectra  and  the  resulting  spectra  are  logarithmically  averaged  to  obtain  the 
initial  estirnale  of  the  neiwu’K-aveiagad  P  wave  spectrum  for  the  current 
event.  The  result  i .  displayed  (-olid  line)  together  with  the  network-aver¬ 
aged  spectrum  expected,  on  average,  for  an  explo-sion  of  that  mj,  value  on 
the  basis  of  [last  experience  at  that  test  site  (dotted  line). 


Compare 


This  module  pernv.!  -  the  operator  !o  furlher  evaluate  the  consis¬ 
tency  of  the  speeira!  data  'iiroiigh  comparisons  of  the  derived  network- 
averageJ  sp  ■etrum  uiou>'d  Imei  with  the  indixulual  station-corrected  spec- 


Reference 


Mi  i,nitu(te  Mtasuremera 


A  -  derived  network-averaged  spectrum 
B  -  individual  station-corrected  spectra 

ira  (solid  lines)  which  went  into  the  network  average.  These  comparisons 
are  displayed  four  stations  to  a  page,  with  subsequent  pages  accessed  using 
the  Continue  function.  During  this  review,  any  stations  showing  significant 
spectral  discrepancies  with  respect  to  the  average  can  be  eliminated  by 
sequentially  selecting  the  Delete  and  station  buttons  from  the  Stations  win¬ 
dow.  When  Continue  is  .selected  following  the  la.st  station  comparison,  the 
Compare  function  is  automatically  terminated.  At  this  point  the  operator 
has  the  option  of  returning  to  the  Avgspec  function  to  recompute  the  net¬ 
work-averaged  spectrum  without  any  stations  which  may  have  been 
deleted  in  the  compare  procc.ss.  In  the  resulting  display,  the  current  net¬ 
work-averaged  specu'um  is  .shown  as  a  .solid  line,  the  previous  estimate  as  a 
dashed  line  and  the  predicted  .spectrurr;  based  on  m^,  again  as  a  dotted  line. 
This  process  of  deleting  (or  adding)  stations  and  recomputing  the  network 
.  average  can  be  repeated  until  a  satisfactory  version  of  the  network-aver¬ 
aged  P  wave  spectrum  is  obtained  for  the  event  under  consideration. 


Reference 


Magnitude  Measurement 


M(Pspec) 

'riiis  module  pcnuiu  ihe  i.pe.  i  ,  s/u.a,,,  a  P  v/ave  spc'.  tral 

magnitude,  M(Pspec),  from  the  computed  u^twvtrk-avcraged  spectrum.  In 
the  accompanying  display,  the  tigure  on  the  left  shows  the  final  comparison 
from  the  Avgspec  module,  while  the  figure  on  the  right  shows  the  compari¬ 
son  between  the  same  observed  specuum  (solid  line)  together  with  the  pre¬ 
dicted  spectrum  (dotted  line)  ciM  ie.sponding  to  the  nij,  value  which  provides 
the  best  overall  tit  to  the  observed  spccU'um  over  the  frequency  band  from 
0.5  to  2.25  Hz.  This  value  of  m^  is  subsequently  designated  M(Pspec). 


Save 

This  tunetion  can  be  ii.scd  by  the  operator  m  replace  the  value  of 
M(Pspec)  which  is  automatically  written  to  the  current  event  magnitude 
file  following  the  first  execution  of  the  M(Psjn\  )  function.  That  is,  if  the 
observed  network-averaged  spectrum  is  subsequently  modified  by  adding 
or  deleting  Ntations,  the  coriesponding  MtPspec)  value  can  be  substituted 
into  the  magnitude  file  using  this  menu  option. 
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'^ok'ciion  ol  this  niodulc  ini'iak's  a  process  by  which  the  operator 
car!  . rTiinJel-hasoil  estimates  of' explosion  yield  and  pP  el'fecLs  corre- 
ii!  :i '  i>)  the  derived  network  avonieed  P  wave  spectrum  for  the  current 
•:  litia!  displiiv  shmvn  in  Pieu'e  23  ol'C'hapter  2  provides  a  menu 
H'le  ti>  n  ut  three  alternate  siiuice  modi  ls:  Mucller/M  trphy  saturated  tuff/ 
:.'  M  Wet  'I'ali'i.  Muei^er/Muirhy  granite  (M/M  Granite)  and  von 
Seeeein'Blandrord  granite  t  VonSee-Bl.mdt ).  Selection  of  any  one  of  these 
itiree  si'urce  models  produces  a  display  showing  a  comparison  between  the 
>'  '!  ii!.'!'  _,'irected  cihs-rved  spe  Ln:m  (solid  line)  and  the  predicted 

speetruin  corresponding  to  the  yield  which  provides  the  best  overall  fit  to 
ihe  !hservati('ns  over  the  band  0 d'O  to  2.25  Hz  (dotted  line).  This  best-fit 
yield  v.ilue  is  listed  at  the  top  oi  the  display  together  with  its  as.sociated 
upper  b('und  uncertainiv  ranee  lection  of  Continue  then  results  in  a  dis¬ 
play  shf'wing  a  compaiisim  between  the  attenuation  and  source  normalized 


'  ii/’per  bound  u/n  rr’iioiry  > 

'  ‘  fit  yield  value 

attenuation-cot  rev  ted  observed  spectrum 

predicted  spectrum  i  onespondim^  to  the  best-fit  \ield  value 


Reference 


Magnitude  Measurement 


A  -  attenuation  and  source  normalized  spectrum 
B  -  predicted  pP  spectrum 


spectrum  (solid  line)  and  the  spectral  modulation  pattern  predicted  by  the 
automatically  determined  best-lit  pP  model  (dotted  line).  The  automatically 
determined  model  parameters  consisting  of  the  ratio  of  pP  to  P  amplitude 
(A)  and  pP-P  delay  lime  (T„),  are  listed  in  parenthesis  in  the  legend.  By 
again  selecting  Continue,  the  operator  initiates  a  process  by  which  this 
inferred  pP  modulation  .spectrum  is  divided  out  of  the  atienuaiion-correcied 
spectrum  to  obtain  the  final  normalised  ob.served  spectrum  which  is  then 
di.splayed  (solid  line)  together  v/iih  the  prediclei)  spec'rum  corre.sponding 
to  the  original  besl  lil  yield  idniicd  line),  flu  inoii;iud  yield  upper  bound 
uncertainty  range  corresponding  to  the  pP-corrected  spectrum  is  also  com¬ 
puted  at  this  time  and  sub.stiluted  into  the  parentiiesis  at  the  top  of  the  dis¬ 
play.  Also  appearing  on  this  display  is  a  box  in  the  upper  right  hand  corner 
with  columns  labeled  A.  T„.  W.  This  indicates  the  option  for  the  operator  to 
interactively  modify  the  automatic  yield  and  pP  estimates.  Thus,  selecting 
the  minus  (plus)  button  under  W  causes  the  yield  to  decrease  (increase)  in 
increments  of  live  percent  liom  the  automatically  determined  yield  and  the 


7« 


A  -  modified  yield  upper  hound  uncertainty  ran}>e 
B  -  interactive  box  used  to  modijy  the  automatk  yield  and  pP  estimates 
C  -  final  normalized  observed  spectrum 

D  -  predicted  spectrum  corresponding  to  the  best-fit  yield  value 

associated  predicted  spectrum  to  be  automatically  updated  in  the  display  to 
reflect  this  new  yield  value.  By  again  selecting  Continue,  the  operator 
causes  a  new  pP  fit  corresponding  to  this  revised  source  estimate  to  be 
determined  and  displayed.  At  this  stage,  the  operator  can  interactively 
modify  the  automatically  determined  pP  parameters  by  selecting  the  plus 
or  minus  buttons  under  the  columns  A  and  T^,  and  the  display  will  continu¬ 
ously  update  to  show  the  comparisons  corresponding  to  these  revised 
parameters.  Once  the  pP  model  has  been  finalized,  the  operator  can  select 
Continue  once  again  to  obtain  the  spectral  comparison  display  correspond¬ 
ing  to  the  interactively  determined  yield  and  pP  parameter  values.  This 
interactive  modification  cycle  can  be  repeated  indefinitely,  thereby  allow¬ 
ing  the  operator  to  explore  a  variety  of  alternative  solutions.  Once  this  pro¬ 
cess  has  been  completed,  the  operator  can  .select  Quit  to  terminate  the 
Yieldfit  function.  At  this  point,  the  operator  can  cither  select  Yieldfit  again  to 
expli  re  a  different  .source  model,  or  .select  Quit  again  to  terminate  the 
Spectral  P  processing. 
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Moment  lensor 


This  module  is  used  lo  eslimatc  wave  momeni  tensoi  solu¬ 

tions  and  associated  explosion  moment  magnitudes  (Mj,)  through  process¬ 
ing  of  the  recorded  long-period  Rayleigh  and  Love  wave  data  from  the 
current  event.  Selection  of  Moment  Tensor  from  the  Magnitude  Measure¬ 
ment  menu  initiates  a  phase  matched  hltering  process  with  associated  dis¬ 
play  and  menu.  Since  it  takes  .some  time  to  initialize  this  process,  it  is  done 
in  the  background  and,  when  completed,  an  icon  appears  on  the  screen 
which  can  be  accessed  by  the  operator  at  his  discretion  to  activate  the  ini¬ 
tial  phase  matched  filter  display  shown  below.  Proceeding  from  top  to  bot¬ 
tom,  this  initial  display  contains  the  menu,  an  information  bar,  a  standard 
analyst  station  display  of  one  component  of  the  surface  wave  recording 
from  the  nearest  station,  the  time  history  output  resulting  from  applying  a 
phase  matched  filter  (Pmf)  to  that  surface  wave  recording  and  finally,  at  the 
bottom  of  the  display  ,  the  amplitude  (left)  and  phase  (right)  spectra  corre- 
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D  -  Pmf  lime  history  output 
E  -  amp! nude  and  phase  spectra 


Reference 


Muiffiitude  Meaxurement 


spending  to  that  Pmf  (Output.  The  Pmf  output  time  history  is  plotted  as  a 
funtuon  of  reduced  time,  computed  with  respect  to  the  arrival  sequence 
predicted  by  the  estimated  group  velocity  dispersion  curve  for  the  propaga¬ 
tion  path  to  that  station.  If  the  data  are  good  and  the  estimated  dispersion 
curve  is  accurate,  this  trace  should  be  dominated  by  a  narrow  pulse  cen¬ 
tered  on  zero  lag  time.  Tlie  amplitude  and  phase  spectra  are  displayed  as  a 
function  of  frequency  over  the  band  extending  from  0.010  to  0.075  Hz  (i.e., 
over  the  period  range  extending  from  about  13  to  100  seconds).  The  aver¬ 
age  value  of  the  amplitude  spectrum  over  this  range  defines  the  surface 
wave  moment  estimate  for  that  station,  which  is  automatically  determined 
by  the  system  and  listed  in  the  infomtation  line.  Ideally,  the  associated 
phase  spectrum  should  oscillate  around  constant  values  of  either  0  or  tt, 
denoting  normal  (sign  +  1)  or  reversed  (sign  -  1)  phase,  respectively.  More 
generally,  in  the  automatic  proces.sing  mode,  the  phase  is  defined  to  be  nor¬ 
mal  if  the  absolute  value  of  the  average  phase  over  this  period  range  is  less 
than  Jt/2  and  defined  to  be  reversed  if  this  quantity  is  greater  than  idl,  and 
the  corre.sponding  value  is  then  listed  in  the  information  line.  As  will  be 
de.scribed  below,  the  operator  can  interactively  re.set  these  automatic  phase 
determinations  in  cases  in  which  the  data  are  ambiguous. 


Pmf 
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\c\t 

Prov  ions 
\ul o  i^il 
Dololo  R«’slotv 
(  Sij;n 


The  functions  listed  under  this  menu  item  permit  the  operator  to 
interact  with  the  displayed  surface  wave  data  and  associated  Pmf  results. 
Selection  of  the  Pmf  option  initiates  an  automatic  process  by  which  the 
moment  and  phase  corresponding  to  the  displayed  waveform  are  deter¬ 
mined  and  written  to  both  the  information  line  and  a  temporary  file.  Selec¬ 
tion  of  Next  (Previous)  refreshes  the  display  to  show  the  Pmf  results 
corresponding  to  the  next  (preceding)  station  in  the  surface  wave  process¬ 
ing  queue.  Similarly,  selection  of  Auto  Pm/ initiates  a  process  by  which  all 
the  usable  surface  wave  data  from  the  current  event  are  automatically  pro¬ 
cessed  through  the  Pmf  routine  and  the  results  are  sequentially  di.splayed 
on  the  screen  without  any  operator  intervention.  Selection  of  Delete/ 
Restore  cau.ses  the  Pmf  results  for  the  displayed  trace  to  be  either  removed 
or  added  to  the  temporary  file  of  pnKessing  results  (depending  on  whether 
it  currently  is  or  is  not  in  that  file)  which  is  to  be  used  as  input  to  the 
moment  tensor  inver.sion  .system  (MTIS)  module.  Finally,  selection  of  the 
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Magnitude  Measurement 


Change  Sign  option  cauMi;.  i  1  1^  .  l.C  JiSplaVCd  tjacc  to 

be  reversed  (i.e.,  to  either  ±  1,  depcndiiie  or  tiic  euneiU  setting). 


MTIS 


Selection  of  this  option  initiates  a  secondarv  process  with  associ¬ 
ated  menu  which  permits  the  operator  to  interactively  process  the  individ¬ 
ual  station  moments  and  phases  determined  in  the  Pmf  module.  The  initial 
display  shows  the  fit  to  the  Rayleigh  (left)  and  Love  (right)  wave  PMF 
moment  estimates  corresponding  to  the  automatic  moment  tensor  solution 
listed  in  the  information  line  below  the  menu  The  default  solution  for  the 
tectonic  component  corresponds  to  thrust  motion  on  a  fault  dipping  at  45° 
and  the  resulting  moment  tensor  solution  is  characterized  by  an  isotropic 
moment  (Ml),  together  with  inferred  strike  and  F  factor  estimates  for  the 
tectonic  component.  On  these  displays,  open  and  closed  symbols  denote 
normal  and  reversed  phase  observations,  respectively.  The  associated  menu 
items  on  this  MTIS  display  permit  the  operator  to  interact  with  these  data 
to  eliminate  questionable  points  and  explore  alternate  solutions. 


A  -  MTIS  menu  hut  tons 
B  -  information  bar 
C  -  nwment  tensor  solution 
D  -  single  station  observations 


Refer cv.'  c 


Mat;nili4de  Measurement 


Eoii 


)  his  function  permits  the  operator  to  select  points  with  the  mouse 
(causing  the  points  to  be  highlighted)  and  then  to  either  change  their  sign  or 
delete  them  from  the  inversion  input  file  by  selecting  the  Change  Sign  or 
Drf  tf  Points  menu  options,  respectively. 


View 


This  function  permits  the  operator  to  either  display  the  seismo- 
gtj'iis  corresponding  to  .stations  (points)  previously  selected  with  the 
mou.se  (Seismogram)  or  to  create  a  graphical  representation  of  the  currently 
displayed  moment  tensor  solution  and  add  it  to  the  Overlays  menu  of  the 
.Satellite  Image  module  for  display  on  the  SPOT  image  of  the  test  site 
{SPOT  image). 

Options 


This  function  permits  the  operator  to  undo  any  previous  selection  of 
points  (Unselect),  to  modify  the  dip  and  rake  of  the  tectonic  solution  from 
the  default  values  for  45°  thrust  (Set  Dip/Rake)  or  to  recalculate  the 
moment  tensor  solution  to  account  for  any  operator  specified  changes  in 
sign  or  deletions  of  points  (Recalculate). 

Selection  of  Set  Dip/Rake  activates  a  display  of  sliders  which  per¬ 
mit  the  operator  to  interactively  set  the  values  of  these  variables  to  any 
phy.sically  permissible  values.  Subsequent  selection  of  Recalculate  from 
the  Options  menu  initiates  a  di.splay  of  the  moment  tensor  solution  corre- 


d,  sliders  used  to  adjust  dip  and  rake  values 
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MCf,.ttiude  Measurement 


Spending  to  these  revised  csiiniaiej  oi  the  dip  ^iiJ  ..i  i',e  niierred  tec¬ 
tonic  solution.  When  the  operator  is  saiisiicd  ihai  a  stable  esiin,  ■  <*  of  the 
moment  tensor  soludon  has  been  obtained,  i(  eari  . k  .d  im  il  ;  ei  tent 
event  database  together  with  the  associated  moment  magnitude  (ivi,j;  using 
the  MTIS  Save  function. 


Save 


Selection  of  this  Pmf  option  inidates  a  process  by  which  the  single 
station  moment  and  phase  estimates  obtained  during  the  current  processing 
session  are  saved  to  the  database  for  subsequent  inpu’  to  the  MTIS  moment 
tensor  inversion  module. 


X4 


Yield  Estimation 


^^1^  Referent 


Yif^ld  Estimation 

This  module  enables  the  operator  to  obtain  unified  seismic  yield 
estimates  and  associated  uncertainty  bounds  for  the  current  event  by  apply¬ 
ing  a  variety  of  statistical  models  to  the  corresponding  network-averaged 
magnitudes.  At  least  one  network-averaged  magnitude  must  have  already 
been  estimated  before  this  module  can  be  activated.  Selection  of  this  option 
from  the  main  FUNCTION  menu  proiluces  a  display  listing  six  different 
stat'stteal  models  and  providing  slider  bars  to  adjust  the  relative  weights  of 
the  ;  iui  magnitude  measures  which  currently  can  be  used  in  the  unified 
yield  analysis  (i.e.,  mb,  M(Lg),  Mo,  M(Pspec)).  These  weights  must  sum  to 
1.00  and  the  initial  default  values  provide  for  equal  weighting  of  all  avail¬ 
able  magnitudes. 


A 


B 


A  -  statistical  models 
R  -  sliders  used  to  adjust  the  w  eifthts 

The  first  five  of  these  models  provide  unified  yield  estimates  and 
extremal  confidence  limits  on  those  estimates  corresponding  to  different 
^  "s  of  constraints  on  a  .specified  nij/yi^dlW)  relation  which,  for  explo¬ 
sions  at  the  Shagan  River  and  Novaya  Zemlya  test  sites  has  been  taken  to 
be 
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Yield  Estimation 


ittfj  =  4.45  +  0.75  log  W 


with  associated  parameters  mj,  (100),  b  and  max  p  where  m5(100)  denotes 
the  estimated  uncertainty  in  the  average  value  of  ni|,  for  explosions  with 
yields  of  100  kt,  b  denotes  the  estimated  uncertainly  in  the  slope  of  the  lin¬ 
ear  relation  between  mj,  and  log  W  and  max  p  denote.',  die  upper  bound  on 
the  absolute  value  of  the  correlation  of  the  yield  estimation  errors  associ¬ 
ated  with  the  different  magnitude  measures.  The  five  available  models  of 
this  type  provide  estimates  conesponding  to  different  values  of  these  three 
parameters  which  are  consistent  with  the  current  range  of  expert  opinion. 
Model  number  6,  labeled  AFTAC,  provides  a  unified  yield  estimate  based 
on  m^  and  M(Lg)  alone,  using  the  current  AFTAC  standard  procedure. 
When  this  model  is  selected,  the  weights  automatically  adjust  to  0.5  for  m^ 
and  M(Lg)  and  to  zero  for  and  M(Pspec). 

Selection  of  any  one  of  these  models  generates  a  display  showing 
the  yield  estimates  corresponding  to  the  individual  magnitudes  plotted  with 
respect  to  the  unified  yield  estimate  and  its  associated  uncertainty  bounds 
(horizontal  lines).  Yield  estimates  determined  from  magnitude  values  not 
included  in  the  unified  yield  calculation  are  shown  as  open  circles  for  pur- 


A 

B 

C 

D 


A,C  -  uncerrainty  hounds 
B  -  unified  }  I  eld 

D  -  value  not  included  in  yield  calculation 
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poses  of  comparison.  The  unified  yield  estimate  (W)  and  its  associated 
uncertainty  factor  (F)  are  listed  in  the  upper  right  hand  comer  of  this  dis¬ 
play  together  with  the  yield  range  corresponding  to  this  F  value  (in  paren¬ 
thesis).  Subsequent  selection  of  models  from  the  menu  produce 
corresponding  displays  beneath  that  for  the  first  choice,  which  remains  in 
place  throughout  the  processing  session  (See  Figure  27  in  Chapter  2). 
When  a  satisfactory  unified  yield  estimate  has  been  obtained,  it  can  be 
saved  to  the  current  event  file  using  the  Save  option  which  appears  on  the 
menu  at  the  top  of  each  display. 
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Statistical  Summary 


This  module  provides  aecc5>s  lo  ii  Ol  stalislical  tests  ano  J;5- 

plays  which  can  be  used  by  the  operator  to  evaluate  the  yield  estimate  for 
the  current  event  in  the  context  of  past  experience  at  that  test  site  and  to 
assess  whether  this  estimate  is  consis'ent  with  the  150  kt  threshold  of  the 
Threshold  Test  Ban  Treaty  (TTBT).  A  unitied  yield  estimate  must  have 
already  been  determined  for  the  current  event  beforts  this  module  can  be 
activated.  Selection  of  this  option  from  the  main  FUNCTION  menu  pro¬ 
duces  a  menu  display  listing  the  available  choices  of  statistical  displays  and 
tests,  and  providing  two  sets  of  slider  bars  which  can  be  used  to  vary  the 
m^/yield  relations  used  in  the  computations.  The  slider  bars  at  the  lower 
left  of  the  display  determine  the  mjj/yield  relation  to  be  used  in  all  tests 
except  the  second  Mixture  of  Normals  test  and  permit  the  interactive  deter¬ 
mination  of  the  intercept  (AO),  slope  (BO)  and  intercept  uncertainty  (AO 
sigma).  Initial  default  values  for  these  parameters  are  4.45,  0.75  and  0.06, 
respectively,  consistent  with  the  nominal  mj/yield  relation  for  explosions  at 
the  Shagan  River  and  Novaya  2iemlya  test  siie.s.  The  slider  bars  at  the  lower 
right  of  the  display  determine  the  mt/yield  relation  for  the  second  Mixture 
of  Normals  calculation  and  have  default  values  consistent  with  the  nominal 
mj/yield  relation  for  explosions  below  the  water  table  at  the  Nevada  Test 
Site  (i.e.,  A1  =  3.94,  B1  =  0.81  and  A1  sigma  =  0.06). 


A  -  statistical  tests 


B  -  slider  bars  used  in  set  the  mj/yiidd  relation 
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Histograms 

Selection  of  either  of  the  first  two  menu  items,  Histogram  of  Histor¬ 
ical  Yields  and  Historical  Yields  versus  Time,  results  in  displays  showing 
the  unified  yield  estimate  for  the  current  event  relative  to  previous  experi¬ 
ence  at  that  test  site,  as  a  function  of  yield  and  dale,  respectively. 

Single  Event  Tests 

The  various  single  event  tests  provide  access  to  formal  compliance 
assessment  tests  which  can  be  used  to  make  quantitative  statements  about 
the  probability  that  the  current  event  is  or  is  not  seismically  compliant  with 
the  150  kt  threshold  level  of  the  TTBT.  Single  Event  Test  1(A0,B0)  assumes 
that  only  indirect,  seismic  estimates  of  the  mj/yield  relation  are  available 
and  utilizes  the  specified  values  of  AO  and  BO  in  estimating  the  compli¬ 
ance  probability.  For  Shagan  River  explosions.  Single  Event  Test  2-JVE 
can  be  acce.ssed  to  factor  in  the  calibration  results  from  the  JVE  CORR- 
TEX  yield  determination  and  utilizes  a  modified  AO  value  which  accounts 
for  this  additional  information  in  the  test  of  compliance.  Similarly,  Single 
Event  Test  3-Exchange  Events  for  Shagan  River  explosions  incorporates 
the  announced  yields  for  the  JVE  exchange  events  into  the  estimation  of 
the  AO  value  used  in  the  compliance  assessment.  The  Single  Event  Test  - 
AFTAC,  on  the  other  hand,  employs  the  current  AFTAC  mt/yield  relation 
and  associated  unified  yield  e.stimate  to  assess  seismic  compliance  of  the 
current  event  with  the  TTBT  ihrc.shold  based  on  just  the  mj,  and  M(Lg) 
magnitude  measures. 

Selection  of  any  of  these  single  event  tests  results  in  a  display 
showing  the  disU'ibution  function  associated  with  the  current  event  unified 
yield  e.stimate  in  relation  to  the  150  ki  threshold.  In  this  representation,  the 
number  at  the  peak  of  the  distribution  corresponds  to  the  unified  yield  esti¬ 
mate  and  the  red  shaded  area  under  the  curve  provides  a  measure  of  the 
probability  of  a  violation  of  the  TTBT  limit.  The  formal  compliance  assess¬ 
ment  resulLs  are  presented  in  tabular  form  in  the  upper  left  hand  comer  of 
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A  -  compliance  assessment  results 
B  -  unified  yield  estinmte 


this  display.  Here,  the  uppermost  box  provides  a  summary  of  the  results  of 
the  hypothe.sis  test  and  indicates  whether  the  current  event  yield  is  or  is  not 
greater  (at  the  95%  confidence  level)  than  the  150  kt  threshold  at  the  nomi¬ 
nal  false  alarm  rale  of  2.5%.  The  first  line  below  that  box  indicates  what 
false  alarm  rate  would  have  to  be  accepted  to  be  able  to  conclude  that  the 
true  yield  of  the  current  event  is  greater  than  150  kt.  The  lines  below  that 
list  the  unified  magnitude  (M)  and  associated  uncertainty  (Osejs).  and  the 
intercept  (AO)  and  slope  (BO)  of  the  mi/yicld  curve  used  in  the  compli¬ 
ance  calculation. 

Mixture  of  Normals 


The  Mixture  of  Nonrials  te.sts  {Mixture  of  N(>rnuils(AO.BO)  and  Mix¬ 
ture  ofNormals(Al,Bl))  provide  more  formal  means  of  comparing  the  uni¬ 
fied  yield  estimate  of  the  current  event  with  the  di.siribution  of  yields  of 
previous  explosions  at  that  test  .site.  In  this  case,  it  is  assumed  that  the 
yields  of  previous  explosions  are  clustered  around  a  few  discrete  values 
and  their  joint  distribution  is  represented  as  a  superposition  of  normal  dis¬ 
tributions,  the  number  and  means  of  which  arc  estimated  by  applying  for- 
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mal  siatistical  criteria.  This  test  can  be  run  using  either  the  (AO.BO)  or 
(  A  I  B! !  definitions  of  the  mj/yield  relation.  In  either  case,  the  resulting 
display  shows  the  unified  yield  of  the  current  event  (*)  with  respect  to  the 
I  .SO  kt  threshold  and  the  best-fitting  mixture  of  normals  solution  for  the  his¬ 
torical  data,  where  the  yields  of  the  individual  distributions  in  the  solution 
arc  listed  beside  their  respective  peaks.  The  tabular  data  in  the  upper  left 
hand  corner  of  the  display  li.sts  the  m|^>  ield  parameters,  number  of  histori¬ 
cal  events  and  individual  di.stribution  a  values  used  in  the  mixture  of  nor¬ 
mals  calculation.  If  both  Mixture  of  Normals  tests  have  been  run,  .selection 
of  the  Mixture  of  Normals  Comparison  option  results  in  a  similar  display  in 
which  color-coded  representations  of  the  mixture  of  normals  solutions  cor¬ 
responding  to  both  the  (AO,BO)  and  (A1,B1)  mj/yield  relations  are  shown 
in  relation  to  the  current  event  yield  estimate  and  the  150  kt  threshold  of 
the  TTBT. 
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Spreadsheet 


This  module  provides  a  less  torinal  environmeni  for  evaluating  the 
seismic  yield  estimate  for  the  current  event  in  the  context  of  the  results 
obtained  from  analyses  of  previous  explosions  at  that  test  site,  as  well  as 
any  available  independent  yield  calibration  data  (e.g.,  CORRTEX).  Selec¬ 
tion  of  this  option  from  the  main  FUNCTION  menu  produces  a  standard 
spreadsheet  display  listing  the  four  individual  seismic  yield  estimates 
(based  on  mj,,  M(Lg),  M^,  M(Pspec))  and  corresponding  unified  yields 
(W(UniO)  for  selected  explosions  at  the  designated  test  site.  For  explosions 
at  the  Shagan  River  test  site,  available  calibration  yields,  including  the  JVE 
CORRTEX  value,  the  JVE  exchange  yields  and  the  yields  published  in  the 
open  literature  are  listed  in  the  column  labeled  W(Bocharov).  The  corre¬ 
sponding  ratios  of  unified  seismic  yield  estimates  to  calibration  yield  val¬ 
ues  are  listed  in  the  column  labeled  Wu/We  and  the  logarithmic  average 
value  of  this  ratio  is  listed  at  the  bottom  of  the  column  as  Avg  Bias, 
together  with  the  estimate  of  the  95%  bound  on  the  observed  ratios  (F). 
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Spreadsheet 


Operator  intciaeiion  with  this  spreadsheet  display  is  restricted  to  the  top 
three  r  'ws  of  entries  specifying  the  magnitude/yield  relations  (Slope,  Inter¬ 
cept)  lor  each  magnitude  type  and  tlic  relative  weight  (Weight)  to  be 
applied  to  each  magnitude  in  the  determination  of  the  unified  yield.  Default 
values  for  the  slopes  and  intercepts  are  4.45  and  0.75,  respectively,  for  m^, 
M(Lg)  and  M(Pspec)  and  5.00  and  1.00  for  Mq.  The  default  values  for  the 
weights  call  for  equal  weighting  of  the  four  magnitudes  and,  however  they 
are  sub.sequently  changed  by  the  operator,  it  is  necessary  that  they  sum  to 
1.00.  Slopes,  intercepts  and  weights  can  be  changed  by  simply  moving  the 
cursor  to  the  appropriate  cell  using  the  cursor  (arrow)  keys  and  typing  an 
equal  sign  followed  by  the  new  numerical  value  and  then  depressing  the 
Return  key.  This  will  cause  the  spreadsheet  display  to  automatically  update 
to  show  the  results  corresponding  to  the  revised  magnitude/yield  relation. 
When  the  analysis  session  is  completed,  this  function  may  be  terminated  by 
sequentially  typing  the  letters  q  and  n. 
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Event  Summary 


This  module  provides  a  one  page  summary  of  the  results  obtained 
by  applying  the  YES  to  the  digital  seismic  data  recorded  from  the  selected 
explosion.  Listed  results  include  the  date,  origin  tinie  and  location  of  the 
explosion,  the  derived  unified  seismic  yield  and  its  as  'C’ated  uncertainty, 
the  results  of  the  statistical  test  of  compliance  of  the  current  event  with  the 
150  kt  limit  of  the  ITBT,  a  tectonic  release  characterization  of  the  event 
based  on  the  surface  wave  moment  tensor  analysis  and,  if  available,  the 
CORRTEX  yield  estimate  for  the  explosion.  In  addition,  each  of  the  esti¬ 
mated  seismic  magnitudes  is  listed  together  with  its  associated  yield  esti¬ 
mate  and  assigned  weight  in  the  determination  of  the  unified  yield  estimate. 


Dependencies  between 
YES  modules 


The  following  chart  shows  the  main  modules  of  YES  and  their 
interdependencies. 
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